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INTRODUCTION 


Many proposed V/STOL aircraft require propulsion systems that deflect 
the thrust during takeoff, landing, and hovering. This deflection can be accom- 
plished by use of flow turning nozzles downstream of the fans, These nozzles 
must operate efficiently and effectively over wide ranges of internal flow condi- 
tions. An important tool in the design of an efficient nozzle is the capability to 
theoretically analyze the nozzle flow. 

Many of the proposed V/STOL nozzles are three-dimensional (fig. 1). 

Since the capability for a full three-dimensional nozzle analysis does not exist 
at this time, simpler but approximate methods for analyzing these nozzles were 
developed. Even though these nozzles are three-dimensional overall, there are 
sections of the nozzles for which a two-dimensional analysis is an informative 
approximation, The analytical procedure described herein is a two-dimensional 
analysis, 

The procedure developed utilizes three computer programs and is similar 
to the procedure for calculating axisymmetric inlet flows (ref. 1) and two- 
dimensional inlet flows (ref. 2). The chief program is the higher order Douglas 
two-dimensional potential flow program (ref. 3) called 24Y at the Lewis Research 
Center, which calculates the incompressible potential flow about arbitrary two- 
dimensional bodies, The other programs, original at Lewis, are called SCIRCL 
and NOZZLEC, Program SCIRCL generates geometric inputs for 24Y from var- 
ious specified analytical shapes or sets of coordinate points for the nozzle com- 
ponents, Program NOZZLEC takes the basic solutions output by 24Y and com- 
bines them into solutions of interest, and, if desired, applies a compressibility 
correction. Figure 2 1s a schematic representation of these solution steps. 

This paper consists of a statement of the problem to be solved, description 
of each program, sample results, and sufficient documentation, including a test 
case, to enable the user to run the programs, 


SYMBOLS 


A, B, C, D combination coefficients 


M Mach number 
V velocity 
W mass flow 


a incidence angle 


Subscripts: 


Cc control station 

com compressible 

i incompressible 

Sl passage between body 1 and body 2 
52 passage between body 2 and body 3 
00 free stream value 


PROBLEM AND SOLUTION DESCRIPTION 
Statement of the Problem 


Geometric representation. - A three-dimensional nozzle can be approxi- 
mated two-dimensionally by using the profiles at the plane of symmetry of the 
three-dimensional nozzle (fig. 3). The basic assumption for two-dimenstonality 
is that all derivatives in the z-direction vanish, that 1s, 0( )/dz = 0. 

Care should be used when attempting to use a centerbody because in two 
dimensions the centerbody is actually a flow splitter. In most cases, adequate 
results can be obtained by omitting the centerbody. 

Nozzles. - The basic problem to be solved 1s to calculate either the com- 
pressible or incompressible flow in an arbitrary two-dimensional nozzle at any 
combination nozzle mass flow rates W. Wo, or Woo. Figure 4 shows three 
possible combinations of nozzle gepiietnigs and weight flow specifications to be 
analyzed. Figure 4(a) shows the single passage nozzle where weight flow rate, 
Ww, Is Specified. Figure 4 (b) shows a nozzle with dual passages where the two 
passage weight flow rates Ws, and Wes are specified. Figure 4(c) shows the 
Same nozzle as figure 4(b). However, in this case the total nozzle and one pas-— 
sage weight flow rate, W, and We, are specified. 

Inlets and nozzles. - While most of the nozzles problems are internal flow 
problems, occasionally external flows about inlets and nozzles are desired 
(fig. 5). If this is the case the programs can accomplish this by inputting values 
of free stream velocity V,, and nacelle incidence angle, a, along with the re- 
quired weight flows described previously. A detailed description of inlet ealeu- 
lation methods is given in references 1 and 2. 

The two-dimensional flow problem 1s solved in several steps (programs). 

1. Geometric representation (Program SCIRCL) 


2, Incompressible potential flow ba~ 5; uttons (Program 24Y) 
3. Combined solutions with compres _ility correction (Program 


NOZZLEC). 
Each step and its corresponding progran ' scribed in the followtng section. 
Description of Solutic . and Programs 

Geometric representation - Progra <L - The nozzle 1s assumed to 
be two dimensional, The profiles are hb into segments at convenient tan- 
gential points (fig. 6). These profiles are } -rred to as bodies in the program 
nomenclature, Each segment may be deti7 : »y an analytic expression or by a 
set of points. The nozzle walls and outer aces (if any) must be extended far 
downstream (fig. 6) to facilitate obtain: :rate potentral flow solutions in 


the nozzle in the region of interest. Th -..i:etry program SCIRCL prepares 
coordinate-point input for efficient use co. vic potential flow program and also 
prints out information such as curvature, « ' angles, etc., which is useful in 
preliminary screening of proposed body s. 


In addition to the surface points, set« ‘ints in the vertical plane, called 
rakes, are needed at axial locations whe -tocity profiles or streamlines are 
desired. At least one set of rake points  .ccuired, corresponding to the axial 


location where the weight flow 1s specified, as discussed with regard to figure 4. 
Such rakes are called control stations. 

Program SCIRCL generates the coo uinates of the rake points for 24Y, 
Program SCIRCL will also produce a Ca ‘c <> plot of the nozzle geometry and 
rake points. 

Incompressible potential flow basic » visons - Program 24Y. - Program 
24Y 1s the Douglas incompressible pote: at flow computer program for single 
or multiple two-dimensional bodies. Br {!\, * e program utilizes a distribu- 


tion of sources and sinks of initially unk <. strength to represent the body pro- 
file. The continuous distribution 1s app. > .ted by representing the body as a 
finite number of discrete flat elements }b — ‘onstant source strength and char- 
acterized by the mid point of the elemen. —ucecl the control point), This approx- 


imation results in a set of parabolic alge «vc equations for source strength at the 
control points that are solved by matrix me. ods. Velocities at the control points 
and at specified off-body points (rake poini.' ‘re then calculated from the source/ 
sink distribution, 

In the two-passage case, the progran used to obtain five basic solutions 
which are used 1n linear combination in ora; to satisfy the prescribed operating 
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conditions (fig. 4). The first basic solution 1s axial flow, the second is a 90° 
cross flow, the third, fourth, and fifth solutions are vorticity solutions about 
bodies 1, 2, and 3, respectively. In the single passage case, only four solu- 
tions are generated, The first two are the same as the two-passage case while 
the third and fourth are vorticity solutions about the two bodies. 

Combined solution - Program NOZZLEC, - This program combines the 
basic solutions V., j=1, 2, 3, 4, 5 from 24Y into any number of solutions 
of interest. A solution of interest is one having specific values of free stream 
velocity, V_, angle of attack of the free stream, a, and weight flow, Ww, 
through the control stations described previously, Optional inputs for W are 
average axial velocity at the control stations (V_, Vgy> or Vago) or average 
Mach numbers ™_,; Mg); or Mg). If either W or M, are specified, they 
are converted to velocities for use 1n the combination solutions. Temperature 
and pressure must also be specified if other than standard conditions are used. 

The method of combination of the basic solutions is shown in figure 7. A 
compressibility correction, if desired, 1s then applied to the velocities, The 
compressibility correction used is 

p, mae Vo 


Cc 


where the terms on the right-hand side gm obtained from the incompressible 
solution, This correction requires no atégration of the inlet geometry and can 
handle locally sonic or supersonic flows. Flow properties (Mach number, pres- 
sure ratio, etc.) are calculated for eithce the compressible or incompressible 
cases depending on the version desired. The compressibility correction is 
actuated by setting ICOMP1 to 0. The tneompressible version 1s generated by 
setting ICOMPI1 to 1. If two passages ar@ specified, two of the rakes mentioned 
under SCIRCL are used as control stations by NOZZLEC_ A control Station is 
the rake where the average axial velocity of the combined solution 1s specified, 
The possible choices of control stations for two passage cases are shown 1n fig- 
ures 4(b) and (c). If only one passage 1s specified (fig. 4+(a)), then any rake can 
be used as a control station. However, it should be noted that the compressible 
solution 1s most accurate in the vicinity of the control station since the compres- 
sibility correction does not exactly satisfy continuity. 

When the velocity on the surface becomes locally Supersonic, the agree- 
ment between theory and experiment is generally not as good as when the flow 
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remains subsonic, To improve the agreement, an optional supersonic correc- 
tion has been incorporated into the program. The supersonic correction is ac- 
tuated as specifying NX as 1 on input, 


Sample Results 


To illustrate the results of the present calculation procedure, the graphic 
output for the Test Case 1s presented next. 

The graphic output of SCIRCL 1s illustrated in figure 8 and consists of a 
plot of the nozzle geometry and the computational point spacing, and the loca- 
tion of the rakes and their point spacing. In addition to its reference value the 
plot is useful for troubleshooting for mistakes in the geometry input and for 
validating the satisfactory distribution of body points and rake points. 

The graphic output of NOZZLEC is illustrated in figures 9 to 11. Fig- 
ure 9 shows the pressure distribution on the surfaces of body 1. Similar plots 
can be obtained for bodies 2 and 3. 

Figure 10 shows the surface Mach number distribution. The comments 
made above for the pressure plots apply also to Mach number plots. 

Figure 11 shows a plot of the flow field vectors, Whenever rake points 
are specified (see fig. 8) velocity vectors are calculated and can be plotted. 


INPUT AND OUTPUT FORMAT 


SCIRCL Input 


Card Format Column FORTRAN Description 
name 

9 AG 1-54 ARE = title for area plots. 
''CROSS-SEC TIONAL 
AREA, SQ, IN,"' 

9 A6 1-54 EX = title for x-axis, 
''AXTAL POSITION, 
IN. tf 

9 A6 1-54 CURVO = title for curvature 


plots. "CURVATURE 
ON SHROUD"! 
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Card Form at Column FORTRAN Description 
name 
9 AG 1-54 SURFAC = title for x-axis with 


surface distance plots. 
"DISTANCE ALONG 
SURFACE FROM 
DOWNSTREAM END 
OF SHROUD, IN," 


These first four cards above will be unchanged for all runs and can be made 
a part of the execution setup deck, or replaced with data statements. 


Card Format Column FORTRAN Description 
name 
1 9 A6 1-54 TITLE Description of Case 


FOR CALCOMP PLOTTING OF INLET PICTURES 


Card Format Column FORTRAN Description 
name 


2 8F10.2 1-10 xX Length, in plot-inches, 
of X-axis required 


11-20 XMIN Value, in data-inches, 
of far left X-pouint. 

21-30 EXEP Data- inch per plot-inch 
along X-axis. 


31-40 YY Length, in plot-inches, 
of Y-axis required. 


41-50 YMIN Value, in data-inches, 
of bottom Y-point 


o1-60 ORD Data-inch per plot-inch 
along Y-axis (usually 
equal to FXFP), 


Card Format Column 


61-70 


71-80 


3 4X, F10.2, 1X 0-14 
3i1 


16 


17 


18 
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FORTRAN 
name 


ELREF 


AREF 


ALPHER 


IPARA 


[F RST 


ISND 


Description 


The X-values in area 
output data are non- 
dimensionalized by 
ELREF. Default 
value ts 1. 


The areas in area out- 
put data are nondimen- 
stonalized by AREF, 
Default value is 1, 


If only one body 1s in- 
put, ALPHER 1s the 
angle of attack (used 
by 24Y) 


Element geometry flag 
used by 24Y 


= 0 Linear elements 

= 1 Parabolic elements 
First-order terms flag 

= 0 No first-order terms 
= 1 First derivative terms 
= 2 Curvature terms 


= 3 Both first-order 
terms 


Second-order terms flag 
= No second-order terms 


= 1 Second derivative 
terms 


= 2 Curvature squared 


terms 


= 3 Both second-order 
terms 


Card Format 


4 2A6, 4, 
211, 212, 
10X, [1, 
9X, 311, 
212 


(Flag ' A’) 


Flag B 


Flag C 


Flag D 


ALL FLAGS 
are 'on' when 
=1, unless other- 


wise noted (Either 
E or J or neither 


can be on but not 
both) . 


Flag J 


Flag E 


Column 


19 


13-16 
17 


18 


19-20 


21-22 


33 


43 
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FORTRAN 
name 


IFLLL 


IDENT 


PROG 


NO6 
LPNCHO 


IPLOTA 


IPLOTC 


IREAD 


[IAB 


IE DON (1) 


Description 


= (0 No combination so- 
lution will be calculated 
by 24Y 


= 1 A combination solu- 
tion will be calculated 
by 24Y 


Six-character tag for 
case I D, 


24Y (right justified) 


0 


1, Do not save output 
for 24Y on Unit 17. 


Plot area against 
X~position 


-1, Plot curvature ver- 


sus X, +1, Plot curva- 


ture versus S 


0 Obsolete) 


Redo geometry from 


point 
Redo entire geometry 


via direct interpola- 
tion 


Card 


Format 


Flag F 
Flag G 
Flag H 
Flag I 


Skip card 5 if J-0 


ty 
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ARTZ 04 


4F10.2 


I4 


3F8.5, 13, 
4F8.5 


Column 


44 
45 
46-47 
48-49 


1-12 


13-24 


1-4 


1-8 
9-16 


‘) 


FORTRAN 
name 


IREDON (2) 
IREDON(3) 
IREDON (4) 
IREDON (5) 


XAA 


YAA 


XBB 


YBB 


ANBDYS 
DELS 


DELSMX 


XRI 


NRAKE 


XRAK 
YLO 


(Note: There is one card for each rake) 


Description 


LPNCHO for any redo 
IPLOTA tor any redo 
IPLOTC tor any redo 


IRE AD for any redo 


X position of starting 


potnt lor potrtiab redo 


Y postlion of star tiny: 


point for partial redo. 


X position of ending 
point for partial redo. 


Y position of ending 
point for partial redo. 


Number of bodies 


Spacing between points 
in region of interest. 


Maximum spacing far 


from region of interest 


Axial distance at which 
surface distance equals 
zero. 


Number of axial locations 
at which data across the 
passage 1s desired. (Can- 
not be greater than 25). 


Axial location of rake. 


Y value of first point 
(lowest point) on rake 
at XRAK, 


Card Format 
9 4F10.2, 
5F8.2 


Column 


17-24 


25-27 


28-35 


36-43 


44-51 


92-59 
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FORTRAN 
name 


YHI 


NY 


XTRAN 


YTRAN 


XSC ALE 


YSCALE 


TYPBDY 


Description 


Y value of last point 
(highesi point) on rake 
at XRAK, 


Number of points 1a raxe 
at XRAK- Restriction 
2 NY = 200. 
nake points are equally 
spaced, AY, between 
Yu. aac YLO where 
AY - Yori - YLO 

(NY - 1) 


Vaiue o! axial transla- 


tO awe 

Ve.ce o. vertical trans- 
28th re voke 

Va.ue of axial scaling of 


rake 
Value of vertical scaling 
ol rake 


Body number. Howeve-s, 
if there 1S symmetry, 
then any body can be m- 
put as a mirror ir2ge 

of any other body “~" . 
can be accomplishec 2 
setting TYPBDY -- - ~ 
Where Mis the nur «>: 
the body to be crea... > 3 
Nis the number a! -:. 
ocdy to be copted 3 ISB 
is set to the Y velus ot i! _ 
“sre about which body N is 
coe m.rrored No other 


ct 


Card 
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FORTRAN 
name 


Format Column 


11-20 ANSEG 


21-30 DELNEW 


31-40 ANLF 


Description 


input is required for 
this body except for 
ANLF 


= Number of segments 
for the particular body, 
except as stated in 
TYPBDY, 


=-1., Delta S spacing 
is set to original value 
of DELS, 


=0., DeltaSts set to 
value of DELS from 
previous body. 


=+ number, Delta Sis 
set to value of input 
DELNEW. 


= 0 Body 1s a lifting 
body, 1.e., in 24Y a 
vorticity solution about 
this body will be calcu- 
lated. 


= 1, Body is a non- 
lifting body, 1.e., no 
vorticity solution will 
be calculated 


Note: All lifting bodies must be input prior to any non-lifting bodies. 


41-48 XTRAN 


49-56 YTRAN 


07-64 XSC ALE 


Value of axial transla- 
tion of this body. 


Value of vertical trans- 
lation of this body 


Axial scaling factor 
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Card Format Column FORTRAN 
name 
65-72 YSCALE 
73-80 XTMAX 
10 3F10. 2 1-10 ENREED 


Description 


Vertical scaling factor 


Maximum value of X for 
which scaling 1s to be 
applied, 


Code indicating type of 
curve to be fitted 
through given points. 


- Q , for bisuper- 
ellipses. See Table I and 
figure 12 for available 
options. Input1, 2, 3, 
4, 5, or 6 (XIN, YIN) 
points as directed. 


= 1000. Same as = 0 but 
with finer point spacing 
near one end of segment 
(two such segments re- 
quired). Usually used 
to give finer spacing at 
the highlight. The 
superellipse going into 
the highhght and the one 
coming out should have 
this flag 


For bisuperellipses where 
the '1000.' option is to be 

used, the rate at which the 
point spacing, ds, changes 
near one end ds ds 4 - 
(Rate) (ds __,) can be speci- 


fied on input. 


Card 


Format 


Column 


13 
FORTRAN 
name 


Description 


The rate (program name 
= PACE) 1s entered as 
the fractional part of 
ENREED for each seg- 
ment. For example, if 
ENREED were input as 
1000.06, the spacing for 
consecutive points would 
be evaluated as follows: 
DS, = DS,_, - (0.06)DS,_ 
if segment 1s to go from 


1 


large-to-small spacing, 
or: DS, = DS. 4 + 15 
(0.06) DS. 4 if segment 
is going from small-to- 
large spacing. 


If PACE 1s entered as 
zero (i.e., ENREED = 
1000.), the default value, 
0.05, 1s used. 


(PACE < 0.133) 


* The first '1000' super- 
ellipse ON A BODY re- 
duces the point spacing 
as far as possible, down 
to a limit of 2 percent of 
the ds value at the begin- 
ning of the segment. 


* All subsequent '1000' 
superellipses input will 
increase ds as far as 
possible up to the input 
value of DELS. 


Card 


Format 


Column 
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FORTRAN 
name 


Description 


* Any number or types 
of segments may be in- 
put between the first 
and subsequent '1000' 
bisuperellipses, with 
the exception of a nor- 
mal bisuperellipse 
(ENREED=0). 


= 1, 1s a Straight line, 
input 2 coordinates 
(XIN(1), YIN(1), 

XIN (2), YIN(2)) 

(fig. 13(a)). 


The first and last 
straight lines on bod- 
les 2 and 3 and the 
last straight line on 
body 1 will automa- 
tically have their 
spacing increased 
from approximately 
DELS near the re- 
gion of interest to ap- 
proximately DELSMX 
away from the region 
of interest. To get 
this type of spacing 
in the first straight 
line of body 1, 
ENREED must be 
specified as 10 


Card 


Format 


Column 


15 


FORTRAN 
name 


Description 


= 10., Special straight 
line used for initiai 
straight line on lower 
shroud. The straight 
line starts with large 
spacing (DEMSMX) 
and ends with small 
spacing (DELS), 

(fig. 13(b)). 


=-l., fits a lemniscate 
between a straight line 
and a point. Input is 
three coordinates 

(fig. 13(c)). 

= -3., fits a cubic be- 
tween two straight 
lines, Input 4 coor- 
dinates (fig. 13(d)). 


= -4.0, generates a 
segment which is a 
mirrored image of 
all the points from 
(XIN(1), YIN(1)) to 
XIN(2), YIN(2)) about 
the line Y = YIN(3), 
See cards 11 and 12 
for XIN and YIN 
formats. 


= 99., for direct in- 
terpolation option 
over one segment 
(See input instruc- 
tions for card 12) 


Card 


11 


12 


12a 


Format Column 

11-20 
21-30 

6F12.5 1-72 

6F12.5 1-72 

Namelist/ 

$ BODYIN / 

Namelist/ 

$ AUXIN / 
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FORTRAN 
name 


REEDEN (1) 


REEDEN (2) 


XIN (1) 
I=1, 2, 3,6, 4, 5 


YIN (1) 
I=1, 2, 3,6, 4,5 


Z (1) 


DONE 


Description 


(See table I). Input ex- 
ponent of x-term for 
bisuperellipse equation 
Blank for all other seg- 
ment types. 


(See table I). Input ex- 
ponent for y-term of 
bisuperellipse. 


X-coordinate for spe- 
ecified points. 


Y-coordinate for spe- 
ecified points. 

Note. If ENREFD = 99, 
instead of cards 11 and 
12, input the following 
cards. 


Z1is a complex array 
containing the X value 
(in the real part) and Y 
value (imaginary part) 
of each given point along 
the segment. The name- 
list will normally be 
ionger than one card 
The program will use 
the input points to genec- 
ate new points with ..ner 
point Spacing near ru.z.c*s 
ol high curvature 

A logical variabie wr ch 


should be input as 
, PRU 
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Card Format Column FORTRAN Description 
name 
BYPASS = . TRUE, 1f no refine- 


ment of tnput points is 
required. 

Note: If ANSEG=0 and 
TYPBDY #0, skip card 
No. 10, and substitute 
lla for 11 and 12a for 12. 


Input Deck Structure 


Card 
1 
2 
3 
4 
5 only if flag J > 0) 
6 
7 
8 
Number of '8' cards = NRAKES 
Number of '9' cards = ANBDYS 
10 
11 
12 Number of '10-11-12' groups for each '9' card 
= ANSEG 
*If ENREED = 99 on card 10, use lla and 12a 
: instead of 11 and 12 
: *if ANSEG = 0 and TYPBDY #0 on card 9, skip 
10 and substitute 11a and 12a for 11 and 12 
10 
11 
12 Figure 14 shows an SCIRCL input form, reflecting 


the above instructions. 
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SCIRCL Output 


Printed Output 


Input file dump (a list of input cards) 

Case number and title 

Input card 3 

Input card 4 (case I.D. and SCIRCL flags) 

Input card 6 (number of bodies, DELS, DELSMX, and XR 

Total number of points for all bodies should not exceed 500. 

Total for any one segment of a body should not exceed 200. 

An error message will indicate if these limits have been exceeded. 
Total number of off-body points must not exceed 200. 

Total number of rakes must not exceed 25. 


Body 1 segment data, body 2 segment data, and body 3 segment data, 


each segment: 
ENREED (as input) and type of segment 
Data depending on type of segment, 
Straight Lane 
xX X (1) X (2) 
Y Y (1) Y (2) 
Last point data* 


Bisuperellipse 


Exponents 

P as readin X X(1) X (2) X (3) X (6) X (4) 
Q as readin Y Y¥ (1) Y (2) Y (3) Y (6) Y (4) 
P as used A XO 

P as used B YO OMEGA 


where A and B= Semimajor and minor axes of transformed superellipse 


XO and YO = Center of the transformed superellipse 


For 


X (5) 
Y (5) 


OMEGA = The difference (in radians) between the slopes of the end point 
slope lines minus 7/2 (1.e., OMEGA 1s a measure of the non~ 


perpendicularity of the slope lines). 
Number of iterations** 
Iteration data*** 
Last data point* 
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‘Magic triangle' messages refer to the triangle formed by extending the 
superellipses slope lines toward each other and drawing a chordline between 1n- 
put points number (X2, Y2) and (X4, Y4). Input points (X3, Y3) and (X6, Y6) 
must lie within this triangle, or, for certain special cases, may lie in a similar 
triangle on the opposite side of the chordline. 

If input points (X3, Y3) and/or (X6, Y6) fall outside the magnic triangle, the 
following message (or messages) is/are printed: 

This point is below the magic triangle, .."' 

''This point is outside the magic triangle...'' 


Cubic 
x X(1) X (2) X (3) X (4) 
Y Y (1) Y(2) Y(3) Y (4) 
A B C D 


where A, B, C, D are the coefficients of the cubic equation. 
Number of iterations** 

Iteration data*** 

Last point data* 


Lemniscate 
x X (1) X (2) X (3) 
Y YQ) = ¥@Q) ¥ (3) 
THETMX CALC AC ALC 


where THETMX CALC = angle between line 1 - 3 and line 1 - 2. 
ACALC from equation R? = 2(ACALC)* sin 2 6 where 


@= THETMX and R? = [XIN(3) - XIN(2)|* + [yin (3) - YIN(2) 7 


Number of iterations** 
Iteration data*** 
Last point data* 

*Last point data. This is the coordinate point data for the last point of the 
Segment. It is presented here because it is overwritten by the first point of 
the following segment and therefore does not appear in the point-by-point array 
below. 

**Number of iterations is the number of iterations required to achieve a 
Satisfactory point spacing for 24Y and have the calculated last point of the seg- 
ment coincide with the input endpoint (to within a prescribed tolerance of 0.1 
DELS). If this satisfactory spacing and end point matching is not achieved with- 
in 100 iterations, the following message is printed out above the number of 
iterations: 
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''This set of data exceeded 150 iterations. Calculations stopped XBRK 
Y BRK. 
X(1), ¥(1), X(2), ¥(2), X(3), Y¥(3), X(4), Y(4), X(S), Y(5)."' 
***Iteration data 


DELS IN Value of ds at end of previous segment 

DELS Value used to start final iteration for this segment 

DELS OUT Value of ds at end of this segment, to be passed on to the 
next segment 

DSTEST Distance from last calculated point to input segment end- 
point 


FINAL PACE The value of PACE at the conclusion or termination) of the 
point-spacing iteration for bisuperellipses with ENREED 
1000 


"INPUT FOR THE NOZZLE COMBINATION PROGRAM" 


For three body case 


NT(1) = _ , NT(2) = , NT(3) = , NP - 
For two body case 
NT(1) = ,j NT(2)= , NP = 


Body coordinates (a separate set for each body) 
Point number 
X-axial distance 
Y-radial distance 
KAPP A-curvature 
DY /DX-slope 
ALPHA-slope angle in degrees 
S-surface distance measured from frist point of each body 
S-S(2) - surface distance measured from XRI 
DELTAS-distance between points 

Rake Information 

XRAK - axial location of rake 
YLO - Y value of first location on rake 
YHI - Y value of last location on rake 
NDY - number of Y points on the rake, as input 
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File output, - The file output (UNIT 17), written by SCIRCL, is used directly 
as input to 24Y. This file consists of the on-body point coordinates and rake 


points. 


Graphic output. - Standard Calcomp can be produced for each geometry run 


as follows: 


(1) For all cases: an X-Y 'picture' of the nozzle with SCIRCL-generated on- 
body points denoted by the '+' symbol and connected with straight lines; segment 
end points have a large circular symbol; off-body rake points are denoted by a 
small square symbol (see fig. 8). 

(2) For cases with IPLOTC #0: a graph of upper shroud body curvature 
against axial position (PILOTC = -1) or against distance along shroud IPLOTC 


= +1). 

Card Format 

1 3([1, 2X), 1X 
7 A4, 5X, 
9(I1, 2X), 
1X, 11 


Column 


11-38 
44 


24Y INPUT 


FORTRAN 
name 


ID 
ISV 


ILIFT 


TTILE 
IPARA 


Description 


Body number, 


Flag to control the saving 
of geometry data for fu- 
ture use again by the Two- 
Dimensional Program. 


= 0 Do not save data. 


= 1 Save the input geom- 
etry data for future use. 


Lift control flag, ILIFT 
= 1-ANLF (from SCIRCL) 


= 0 This 1s not a lifting 
body. 


= 1 This is a lifting body. 
Body description 
Element geometry flag 

= Q Linear elements 


= ] Parabolic elements 
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Card Format Column FORTRAN 
name 
47 IFIRST 
50 ISEC 
53 ITR 


Description 


First-order terms flag. 
= 0 No first-order terms 
= 1 First derivative term 
= 2 Curvature term 


= 3 Both first-order 
terms 


Second-order terms flag. 


= 0 No second-order 
terms 


= 1 Second derivative 
term 


= 2 Curvature squared 
term 


— 3 Both second-order 
terms 


Geometry transforma- 
tion flag, 


= 0 Transformation card 
will not be input. 


= 1 Geometry transfor- 


mation card will be input, 


= 2 Ellipse generation, 
Ellipse generation card 
will be input. Transfor- 
mation card will not be 
input, 


3 Ellipse generation, 
Ellipse generation card 
will be input. Transfor- 
mation card will be input. 
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Card Form at Column FORTRAN 
name 
56 INORM 
59 IBOD 


Description 


Geometry normaliza- 
tion flag. 


= 0 Geometry will not 
be normalized 


= ] All of the geometry 
data (X and Y) will be 
divided by the chord 
length before use by the 
potential flow program. 


Body disposition flag. 
This flag together with 

the IDOLD parameter 
controls the sequence of 
shapes that are to be pre- 
sented to the potential 
flow analysis part of the 
program, With the useof 
these two flags and the 
ISV parameter it 1S pos- 
sible to perform a variety 
of multi-element analysis 
problems with a minimum 
of input data. For normal 
useage when all the geom- 
etry data are input only the 
IBOD = 1 and =2 inputs are 
used. 


1 New geometry is being 
input. The storage of gco- 
metry data for the polten- 
tial flow solution will start 
with this body. 


Card 


Format 


Column 


62 
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FORTRAN 


name 


IDOLD 


Description 


= 2 New geometry 1s be- 
ing input but this 1s not 
the first body. This body 
will be added to the sc- 
quence of body data al- 
ready input. 


3 New geometry 1s be- 
ing input but it 1s to be 
added to an old sequence 
of data 


= 4 All previously saved 
geometry will be used. 


= 5 The geometry for this 
body will be selected 
from the previously saved 
data (body IDOLD will be 
selected). This selected 
body will be added to the 


current string ol bodies 


= 6 Previously saved geo- 
metry data will be used 
with the body number in- 
dicated by the IDOLD pa- 
rameter removed from 
the solution 


Old body [ID number. This 
parameter 1s used in con- 
junction with the [BOD pa- 
rameter in selecting which 
previously saved shape 1s 

to be retrieved as the pre- 
sent body 


Card Format Column FORTRAN 
name 
65 IPVOR 
68 LAST 
72 ITYPE 
2 This card is input only when ITR = 1 or 3. 
7 (F8.0, 1X), 1-8 CHORD 
F8.0, [1 
10-17 XMULT 
19-26 YMULT 
28-35 DX 
37-44 DY 
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Description 


Vorticity distribution 
flag. 


= 0 Use constant vor- 
ticity between body ele- 


ments. 


= 1 Use variable vor- 
ticity distribution 


Last body flag. 


= Q This 1s not the last 
body. After this body 
1s input the program 
will return to read 
another Body Title and 
Control Card for the 
next body. 


= 1 This 1s the last body. 
=] 


Body chord length. [f 
input as = 0.0 the pro- 
gram will calculate 
CHORD. This parameter 
used only when INORM 

= 1. 


X-coordinate multiplier. 
Y-coordinate multiplier. 


Delta-X increment ap- 
plied to the X-coordinates. 


Delta-Y increment ap- 
plied to the Y-coordinates., 
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Card Format Column FORTRAN Description 
name 
46-53 THETA Angle of rotation (posi- 


tive for a counterclockh- 


wise rotation). 


55-62 XTO X-coordinate of the 
center of rotation. 


64-71 YTO Y-coordinate of the 
center of rotation. 


72 ITYPE - 2 


Note: These transformations are applied in the following order: 
1. Rotate by angle THETA. 
2. Shift by DX and DY. 
3. Apply multiplier factors, XMULT/CHORD, and YMULT/CHORD. 


3 6F10.0, 4X, 1-10 X(1) X-coordinate of the 
I1, 2X, I, 11-20 X (2) geometry. Up to six 
3X, [1 points may be input 


21-30 X(3) on each card depend- 


31-40 X (4) ing upon how the INO 
41-50 X(5) flag is set 

21-60 X (6) 

65 INO Number of data points 


per card, If there are 
6 values per card, INO 
may be left blank. 


68 IST AT Last card flag. 
Q This 1s not the last 
X-coordinate card, 
More cards will lollow, 


1 This is the last 
X-cooridinate card. 


~] 
lw 


ITYPE 5 


Card Format 
4 6F10.0,4X 
I1, 2X,11, 
3X, I1 
5 2X, 13, 5X 
F10.0, 51X, 
Il 
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Column FORTRAN 
name 


1-10 Y (1) 
11-20 Y (2) 
21-30 Y(3) 
31-40 Y (4) 
41-50 Y(5) 
51-60 Y (6) 


65 INO 
68 ISTAT 
72 ITYPE 


Description 


Y-coordinate of the 
geometry. Up to six 
points may be input on 
each card depending 
upon how the INO flag 
is set. 


Number of data points 
per card. If there are 
six values per card, 

INO may be left blank, 


Last card flag. 

= 0 This is not the last 
Y-coordinate card, 
More cards will follow. 


= j This is the last 
Y-coordinate card, 


= 4 


This card is input only when ITR = 2 or = 3. 


3=5 LX 


11~20 ELPSTH 
72 [TYPE 


Number of points to 
be generated, 


Ellipse thickness ratio. 


=95 


Note: At this point in the input the geometry for one body has been read in and 
stored. If another body is to be loaded (i.e., a multiple element case) 
and the LAST parameter on the Body Title and Control Card was = 0, 
then another Body Titl e and Control Card will be expected next. If 
LAST = 1 then the geometry load has been completed and the program 
execution cards will be expected next. 
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Two-Dimensional Program Execution Cards 


All of the following cards are associated with the potential flow solution and 
are input after the geometry data has been loaded. For most normal problems 
only the Flow Title Card and the Flow Control Card are required to complete 
the analysis. All of the other flow execution cards are associated with other 
options normally not used such as non-uniform flow. 


Card Format Column FORTRAN Description 
name 
6 15A4, 11X, 1-60 FTITLE Title or description of 
[I] case 
Le ITYPE & 
7 11, 4X, F10.5, 1 INCLT Cy - @ flag 
2K, 11, 2%, 0 Angle of attack, a, 
P10 5, 5(44, 1S inpul. 
I1), 9X, I1, 4X, 
11, [2, 12 - 1 Total lift coefficient, 
Cr 1S input. 
6-15 CLT Value of angle of attack 


or lift coefficient depend- 
ing upon how the INCLT 
flag was set. 


18 ICHORD Reference length flag, 


= 0 The reference length 
used in calculating the 
Cr is set- 1 0 


1 The reference length 
used in calculating the 
Ce will be input as the 
CCL parameter 


21-30 CCL The input value for the 
reference length (chord) 
used tn calculating the 


Cy, 
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Card Format Column FORTRAN 
name 
35 IND 
40 ISOL 
45 IOFF 
50 NONU 


Description 


Print control flag. 


= 0 Minimum print of 
output solutions 


= 1 Pring the individual 
solutions (0°, 90°, ie 
etc.). 


Matrix solution method 
control flag, 


= 0 Use routine SOLVIT 
for the matrix solution 
(used when a very large 
number of geometry 
points have been input). 


= 1 Use routine QUASI for 
the matrix solution. 


= 2 Use routine MIS1 for 
the matrix solution. 
Maximum number of geo- 
metry points is = 101. If 
number of points is great- 
er than 101 the program 
will automatically shift 

to use SOLVIT. 


Off-body calculation flag, 


= 0 Off-body points will 
not be calculated, 


= 1 Off-body points will 
be calculated. 


Non-uniform flow flag. 


= 0 Non-uniform flow is 
not input. 


Card Format Column 


20 


65 


70 


30) 


FORTRAN 
name 


NBNU 


IPRINT 


MORE 


[Ir LUL 


Description 


#0 Non-un:tform flow 
will be input. The 
number of flows input 
is = NONU (maximum 
of 6 permitted). When 
this option 1s used the 
program automatically 
sets the parameter 
ISOL = 1. 


The number of bodics 


for which the non- 
uniform flows arc input. 


Print/punch flag. 
- 0 Normal output. 


= 2 Print the individual 


matrices. 


7 Punch the output on 
cards. 


Last case flag. 


= 0 This 1s the last solu- 
tion case, 


- 1 This 1s not the last 
solution case Another 
set of Flow Title and 
low Control Cards (and 
any non-uniform or off- 
body cards) will be ex- 
pected next after this 
case ts completed 


( No combination solu- 
lution calculated 


Card Format 
8 6F10.0, 
10X, [2 
9 I1, LX, 301 
10(2X, [3) 
F10.5,5X, 
I2 


Column 


11-20 
etc. 


91-60 


71-72 


ol 
FORTRAN 


name 


ITYPE 
CNU(1) 


CNU(2) 


CNU(6) 


ITYPE 


Description 


= 1 Combination solu- 
tion calculated 


= 9 
Combination constant 
for the first non- 


uniform flow. 


Combination constant 
for the second non- 
untform flow. etc. 


Combination constant 
for the sixth non- 
uniform flow. 


= 10 


This card 1s input only when NONU 7 0. 


IBOD 


IN 


IT 


Body ID number for 
which non-uniform flow 


1S Lnput 


Normal velocity input 
flag. 
= 0 Normal velocities 


are not input, 


= 1] Normal velocities 


are input. 


Tangential velocity 


input flag 


- 0 Tangential veloci- 
ties are not input 


Card 


10 


Format Column 


8-10 

13-15 
18~20 
23-25 
28-30 
33-35 
38-40 
43-45 
48-50 
93-90 
06-65 


71-72 


32 


FORTRAN 
name 


NN 


11(1) 
IE (1) 
11 (2) 
IE (2) 
11 (3) 
LE (3) 
11 (4) 
IE (4) 
11 (5) 
IE(5) 
CB 


ITYPE 


Description 


= 1 Tangential velocities 
are input, 


I] and IE specify the 
element range for which 
the non-uniform flow ts 
input. [1 1s the first 
element of the range 
and IE is the last ele~ 
ment of the range. Up 
to five ranges are per- 
mitted per body. 


Flow scaling coefficient 
for this body. 


ad 


This card 1s input only when NONU # 0 and IN =1. One set of these cards 
must be present for each body for which non-uniform flow 1s to be input 


(when NONU # 0 and IN = 1), 


The number of data points in each set 1s 


determined by the number of element ranges input and by the number of 


elements in each element range. 


6F10.0, 10X, 1-10 
fa 11-20 
21-30 
etc. 


21-60 


71-72 


This card 1s input only when NONU 40 and IT 


VNUF() 
VNUF (+1) 
VNUF (+2) 


ITYPE 


Input values of non- 
uniform normal veloc- 
itv. Input six valves 


per card, 


i 


1. One set of these cards 


must be present for each body for which non-uniform flow 1s to be input 


Card 


Note: 


Card 


12 
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(when NONU #0 and IT = 1). The number of data points 1n each set 1s 
determined by the number of element ranges input and by the number of 


elements in each element range. 


Format Column 


6F10.0, 10X, 1-10 
[2 11-20 


21-30 
etc, 


51-60 
71-72 


FORTRAN 


name 
VTUF() 
VTUF(F1) 
VTUF (+2) 


ITYPE 


Description 


Input values of non- 
uniform tangential 
velocity. Input six 
values per card. 


=13 


The cycle of input for the cards that follow the Flow Control Card 1s con- 


trolled by the NONU and NBNU parameters. 


follows: 


Flow Control Card (NONU and NBNU are input) 


Non-Uniform Flow Coefficients 


Non-Uniform Flow Body Card 
Non-Uniform Normal Velocities 
Non-Uniform Tangential Velocities 


Format Column 


FORTRAN 


name 


The order of cards is as 


these cards 


these cards are | 4re cycled 


cycles for for NONU 
NBNU bodies flows 
Description 


The following cards related to off-body calculations are input only if 


IOFF = 1. 

11, 9X,7A4, 1 
12X, 2(2X, I1), 
2(5X, 11), 2X, 

12 


11-38 


ID 


TITLE 


Identification number for 
this group of off-body 
points. Off-body points 
are read in groups of up 
to 100 maximum at a 
time. There is no limit 
on the number of groups. 
Title or description for 


this group of off-body 
points. 


Card Format Column 


od 


96 


62 


68 


41-72 


34 


FORTRAN 
name 


ITR 


INORM 


IDOLD 


LAST 


[TYP] 


Description 


Coordinate transforma- 
tion flag. See the ITR 
parameter on card 1. 


Coordinate normaliza-— 
tion flag, 


= 0 Off-body coordinates 
will not be normalized 


= 1 Normalize the coor- 
dinates by the input 
chord or by the chord for 
body with ID IDOLD, 


Body selection [lag for 
normalizing olf-body 
points, 


= 0 Use input chord 
(Type 22 Card) to nor- 


malize off-—body points, 


#0 Use chord for body 
with ID = IDOLD to nor- 
malize the olf-body 
points, 


Off-body group termi- 
nation flag. 


Q Additional group(s) 
of off-body points will 
be read in after this 
eroup 1s completed. 


- ] This is the last 
eroup of off-body 
points. 
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Card 


13 


14 


Format 


Column 


30 


FORTRAN 
name 


Description 


This card 1s input 1f ITR = 1 or = 3 on the Off-body Title and Control 


Card. 


7 (F8.0, 1X), 
F8.0, 11 


6F 10.0, 4X, 
Il, 2X, 11, 
3X, [1 


1-8 


10-17 
19-26 
28-35 


37-44 


46-53 


99-62 


64-71 


72 
1-10 
11-20 
21-30 
31-40 
41-50 
01-60 


CHORD 


XMULT 
YMULT 
DX 


DY 


THETA 


XTO 


YTO 


ITYPE 
X(1) 
X (2) 
X (3) 
X(4) 
X (5) 
X(6) 


Body chord length This 
parameter used only when 
INORM = 1 and ID®@LD = 0 


X-coordinate multiplier. 
Y-coordinate multiplier. 


Delta-X increment ap- 
plied to the X-coordinates. 


Delta-Y increment ap- 
plied to the Y-coordinates, 


Angle of rotation (posi- 
tive for a counterclock- 
wise rotation). 


X-coordinate of the center 
of rotation, 


Y-coordinate of the center 
of rotation, 


=2 


X-coordinate of the off- 
body point. Up to six 
points may be input on 
each card, depending 
upon how the INO flag is 
set. 


Card 


15 


Format 


6F10.0, 4X, 
T1, 2X, 11, 
3X, [1 


Column 


65 


68 


72 
1-10 
11-20 
21-30 
01-40 


41-50 
21-60 


65 


68 


ee 
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FORTRAN 


name 


INO 


ISTAT 


ITYPE 
Y (1) 
Y (2) 
Y (3) 
Y (4) 


Y (5) 
Y (6) 


INO 


IST AT 


ITYPKF 


Description 


Number of data points 
per card. If there are 
six values per card, 

INO may be left blank, 


Last card flag 


- Q This 1s not the Last 
X-coordinate card. More 
cards will tollow. 


- 1 This is the last 
X-coordinate Card. 


= 8 
Y-coordinate of the off- 
body point, Up to six 
points may be input on 
each card depending 


upon how the INO flag 1s 


set. 


Number of data points 
per card. [f there are 
six values per card, INO 
may be left blank. 


Last card flag. 


0 This 1s not the last 
Y-coordinate card. More 
cards will follow. 


1 This is the last 
Y-coordinate Card, 


a 


Card 


16 


3? 


Format Column FORTRAN Description 
name 


This card 1s input only when ITR = 2 or = 3 on the Off-Body Title and 
Control Card. 


2X,13,5X, 3-5 LX Number of points to be 
F10.0, 51X, generated. 
ty 11-20 ELPSTH Ellipse thickness ratio. 
72 ITYPE = 
24Y OUTPUT 


Printed Output 


For each body the following data are output: UNTRANSFORMED COORDINATE 
DATA FOR BODY xx and ELEMENT COORDINATE DATA FOR BODY xx. The un- 
transformed data are the geometric points which were input from SCIRCL. For the 
element data, the midpoint of each segment 1s calculated. Also printed are the 
linear distance and surface distance between each segment. For each element, 


the sine, 


cosine, and local curvature are also output. Following the element data 


is the Antal surface distance for that body. 
After the input and transformed coordinates are output a BODY GEOMETRY 
SUMMARY table is output. The following is output in the summary: 


BODY DESCRIPTION Title of test case 

BODY ID Body number 

LIFT TYPE = Yes, vorticity solution about this body is 
calculated 


N/O 
SID 
TFORM 


= No, no vorticity solution calculated 
New or old body data 
= S, body geometry data saved 


= Yes, body coordinates have been trans- 
formed (see ITR 1n input) 


= No, body coordinates have not been trans-~ 
formed 
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NORM = Yes, body coordinates have been normal- 
ized (see INORM 1n input) 


= No, body coordinates have not been normal- 


ized 


CHORD If INORM = 1, the CHORD 1s the body chord 
length. If INORM = 0, CHORD 000 


TYPE Type of solution to be generated. 
= L, only linear terms 
= P, parabolic terms 


SIGMA If TYPE = P, SIGMA describes which para- 


bolic terms considered 


= 1, First derivative term 


F 

C = 1, Curvature term 

S = 1, Second derivative term 
C 


= 1, Curvature squared term 


ELEMENT STORAGE 


N/O See N/O above 
FIRST Number of first element for given body 
NO Number of elements for given body. 


Total Number of Bodies 
Self explanatory 
Total Number of Elements 

Following the GEOMETRY SUMMARY table 1s a table of the vorticity for each 
element. If IPVOR does not equal zero then a variable vorticity for each element 
is calculated If IPVOR equals zero then the vorticity 1s a constant for each ele- 
ment. 

Following the vorticity functions 1s a printout of the combination coefficients 
used internally by 24Y for its version of a combined solution based upon input 
alpha or total lift if INCLT = 1. 

After the combination coefficients of INDIVIDUAL FLOW SOLUTIONS for each 
body are printed. The first solution is axial flow and the second individual solu- 
tion 1s for crossflow. There 1s also an individual solution based upon each body 
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which is a lifting body. Each individual solution contains the following: 


POINT NUMBER Element number 

VN Normal velocity at the point 
VT Tangential velocity at the point 
SIGM A Source strength density 


The individual solutions are also printed on a disc file to be used as input to 
NOZZLEC as basic solutions. 

The next section of output is the COMBINED SOLUTION if IFILL =1. This 
solution is the result of the combination of individual solutions into a solution 
which satisfies the input conditions of either alpha or total lift. Also included are 
the integrated values of lift for each body. 

If there are off-body points specified in the input, the next section is the indi- 
vidual solutions at each off-body point. These solutions correspond to those for 
the on-body points. Also included is a combined solution for the off-body points 
(if IFILL = 1) using the same combination coefficients as those for the on-body 
points. 


NO ZZLEC Input 


English engineering units are used throughout the program. 


Length, in. 
Velocities, ft/sec 
Angles, deg 
Pressure, lb/ft” 
Temper atures, °R 
Densities, slug /ft® 
Force, lb 

Weight flow, lb/sec 


Card Form at Column FORTRAN Description 
name 
1 3.A6 1-18 TITLE Title card 
2 914 1-4 NT Total number of on-body 


points. [NT (3) from 


Card Format 


13-16 


17-20 


25-28 


Column 


40 


rORTRAN 
name 


NS1 


NH 


NP 


NX 


KND 


Description 


SCIRCL (3 body case), 
or NT(2) from SCIRCL 
(2 body case)] . 


Number of on-body 
points on body 1. 
(NT(1) from SCIRCL). 


Number of on-body 
points on bodies 1 and 

2 (NT(2) from SCIRCL) 
If there are only 2 bod- 
tes, NH~- number of on- 
body points on body 1 
only. (NT(1) from 
SCIRCL)., 


Total number of off-body 
points. (NP from 
SCIRCL). 


Flag for type of input 
through the control sta- 
tion: 

0, weight flow 

1, Mach number 


2, velocity 


1, apply supersonic 
velocity correction 
to data 


Flag for scaling variables. 
All tnput lengths are divi- 
ded by ELND 
If KND -1, KLLND YCU 
0, FLND 1 
it, FLND YCU 
- \CL 


Card 


3 


Format 


10F8.0 


Column 


29-32 


33-36 


1-8 


9-16 


17-24 


29-32 
33-40 


41-48 


49-56 
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FORTRAN 
name 


ICOMP1 


IHUB 


VC 


VS1 


V8S2 


VINF 
ALFA 


MC 


MCl 


Description 


= 2, ELND-=- the 
read-1n value 
from card 4, 


Type of solution to be 
computed: 
= 0, compressible 
= 1, incompressible 


IHUB = 0, 2 body case 
#0, 3 body case 


Average axial velocity at 
the control station be- 
tween bodies 1 and 3 up- 
stream of the body 2. 


Average axial velocity at 
the control station be- 
tween bodies 1 and 2. 


Average axial velocity at 
the control station be- 
tween bodies 2 and 3. 
Note 1f IW=2, then two of 
the three velocities must 
be input for the two pas- 
sage case 


Free stream velocity 


Angle between free stream 
velocity and X-axis of the 
nozzle 


Average Mach number at 
the control station be- 
tween bodies 1 and 3, 


Average Mach number at 
the control station be- 
tween bodies 1 and 2. 


Card 


4 


Format 


10F8 .0 


Column 


of-64 


65-72 


73-80 


1-8 


17-24 


tu 

wh 
i 

IN 


ve 


ond 


FORTRAN 
name 


MC2 


TTOTAL 


PT 


ELND 


WDOTC 


WDOTC1 


W DOTC2 


Description 


Average Mach number 
at the control section 
between bodies 2 and 3. 
If IW=1 then two of the 
three Mach numbers 
must be input for the 
two passage case 


Total temperature. If 
TTOTAL = TSTAT=0, 
then TTOT AL = 518.67 
will be used, 


Total pressure, If 
PT=0.0 and PSTAT#0.0, 
the program will calcu- 
late PT, If PT=0.0 
and PSTAT=0.0, PT 1s 
set to 2116. 


ELND 1s the arbitrary 
length used for scaling 
or normalizing. Refer 
to KND input. 


Weight flow at the con- 
trol station between 
bodies 1 and 3. 


Weight flow at the con- 
trol station between 
bodies 1 and 2. 


Weight flow at the con- 
trol station between 
bodies 2 and 3, If IW -0, 
then two ol the three 
weight flows must be in- 
put for the two passage 


case, 


Card Format Column 


33-40 
41-48 


49-56 


57-64 


65-72 


7 3-80 


4A 6F10.0 1-10 


11-20 
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FORTRAN 
name 


PSTAT 
TST AT 


CUTOFI1 


CUTOF3 


CUTO F2 


VPERIN 


XMIN 


Description 


Static pressure 


Static temperature. 

If PSTAT and TSTAT 

are not 0.0, total tem- 
perature and total pres- 
sure will be calculated 
using PSTAT and TSTAT, 


If CUTOF1 7 0, then the 
pressure ratio P./P.7 

on body 1 will be plotted 
against a dimensionless 
surface distance S/CUTOF1 
starting at X = XR1 for a 
distance of S= CUTOF1, 


Same as CUTOF1 except 
for body 3. 


Same as CUTOF1 except 
for body 2. 


If VPERIN 1s greater than 
zero then a CALCOMP 

plot of the nozzle showing 
rake point flow field data 
will be produced. VPERIN 
is the value of a unit vector 
in ft/sec/in. If VPERIN 

is not equal to zero, card 
't4a'' must be input. 


The length in plot inches 
of the abscissa of the v~l- 
ocity plot. 

Value, in data inches, of 
far left X-pouint, 
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Card Format Column FORTRAN Description 
name 
21-30 EXEP Data inch per plot inch 


along X-axis. 


31-40 YY Length, in plot inches, 
of the ordinate. 


41-50 YMIN Value, in data-inches, 
of bottom Y point. 


51-60 ORD Data inch per plot inch 
along Y-axis, 


5 3F10.0 1-10 XTEST Axial location of the 
control station up- 
stream of body 2 be- 
tween bodies 1 and 3. 


11-20 YCL Y on body 1 at XTEST 
21-30 YCU Y on body 3 at XTEST 
6 3F10.0 1-10 XTEST1 Axial location of the 


control station be- 
tween bodies 1 and 2, 


11-20 YCL1 Y on body 1 at XTEST1, 
21-30 YCU1 Y on body 2 at XTESTI1. 
7 3F10.0 1-10 XTEST2 Axial location of the 


control station between 
bodies 2 and 3. 


11-20 XCL2 Y on body 2 at XTEST2 
2136 YCU2 Y on body 3 at XTFST2 


Note: Cards ''6"' and ''7'' are not used for a single-passage case, 


8 3F10.0 1-10 XR1 Axial location on body 1 
where SS 0, 
11-20 XR3 Ax al location on body 33 


where SQ, 


Card Format Column 
21-30 

9 3F10.0 1-10 
11-20 
21-30 
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FORTRAN 
name 


XR2 


YR1 
YR3 
YR2 


Figure 15 1s the input form for NOZZLEC. 


NOZZLEC Output 


Printed Output 


Input file dump 


TITLE - followed by 2-D COMBINATION SOLUTION 
Version of run (i.e., compressible or incompressible) 
A list of the basic flow solutions obtained from 24Y 


Description 


Axial location on body 2 
where S- 0. Leave 
blank for two body case, 


Y on body 1 at XR1. 
Y on body 3 at XR2, 
Y on body 2 at XRH, 


In the table that follows, several functions of four different velocities are given. 


The velocities are: 


Control: VC, average axial velocity at upstream control station 
Lower passage: VS1, average axial velocity at control station between bodies 


1 and 2 


Upper passage: VS2, average axial velocity at control station between bodies 


2 and 3 


Free stream: V,.» free stream velocity 


The rest of the table is self-explanatory perhaps with the exception of the terms 
INC and COMP. INC means calculated from the incompressible equations and 
COMP means calculated from the compressible equations, The rest of the output 


will be defined by its name, 


ALPHA Angle of attack of nozzle 
VINF/VC V/V 3 
VSONIC Critical velocity uncorrected for compressibility 


VSONICC Critical velocity 


WDOTCR 


WDOTLCR 


WDOTUCR 


TSTAT 
PSTAT 
PSTATC 
ASTAT 
RHOSTAT 
WDOTC 
WDOTL 
WDOTU 
VIC 


VICL 


VICU 


TTOT 
PTOT 
PTOTC 
ATOT 
RHOTOT 
THET 
DEL 
XRI1 


YRI1 
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Corrected weight flow at upstream control station 


WDOT xa/ THET 


DEL 


Corrected weight flow at control station between bodies 
1 and 2 


Corrected weight flow at control station bctween bodies 
2 and 3 


Free stream static temperature 

Free stream static incompressible pressure 

Free stream static compressible pressure 

Free stream static speed of sound 

Free stream static density 

Input mass flow at upstream control station 

Input mass flow at downstream lower control station 
Input mass flow at downstream upper control station 


Incompressible average velocity at upstream control 
station 


Incompressible average velocity at downstream lower 


control station 


Incompressible average velocity at downstream upper 
control station 


Free stream total temperature 

Free stream total pressure incompressible 
Free stream total pressure compressible 
Free stream stagnation speed of sound 
Free stream stagnation density 
TTOT/518.67 

PTOTC/2116., 22 

Input 


Input 
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XRI2 Input 

YRI2 Input 

XTEST Input 

YCL Input 

YCU Input 

LND Length used for scaling 

XTEST1 Input 

YCLi1 Input 

YCU1 Input 

XTEST2 Input 

YCL2 Input 

YCU2 Input 

P-S CUTOFFI1 Input 

P - 8S CUTOFF2 Input 

P -S CUTOFFS3 Input 

NT Input 

NP Input 

NS1 Input 

NH Input 

KND Input 

[IW Input 

NX Input 

ICOMP1 Input 

IHUB Input 

V1 Average axial velocity at specified control station for 
basic solution 1 from 24Y 

V2 Same as V1 except for basic 24Y solution 2 

V Same as V1 except for basic 24Y solution 3 


V4 
V5 
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Same as V1 except for basic 24Y solution 4 


Same as V1 except for basic 24Y solution 5 


24Y basic solutions 1 and 2 are axial and crossflow solutions. For a three-body 


case, solutions 3, 4, and 5 are vorticity solutions about bodies 1, 2, and 3, re- 


spectively. For a two-body case, solutions 3 and 4 are the vorticity solutions 


about bodies 1 and 3. 


A, B, C, D 
VINFP 


OTHER MESSAGES: 


Note there is no 5 solution for this case. 
Coefficients of combination 


Incompressible free stream velocity ''uncorrected'' for 
compressibility if the input value was compressible 


"VRESON = IS GREATER THAN VMAX, 
VCONC = .'' The velocity at a certain on- 
body point exceeds the allowable value for the local ex- 


pansion condition so that the isentropic ratio term: 
1. - VCONC 1s less than zero. Where, 


2 roar 


'"'T EXCEEDS 20 ITERATIONS FOR RHOBAR., 

VBAR = » VCOMP = » RHOBAR = 

VBAR HAS BEEN REDUCED TO VCOMP «x 
RHOBAR/RHOTOT"', Subroutine VBARIT attempts to 
find the average density at each axial location using the 
isentropic density ratio, the stagnation density (RHOTOT) 
and the average incompressible VBAR (based on weight 
flow and the cross section). It has failed. VCOMP 1s the 
20th attempt at finding the compressible velocity and has 
been used to compute the RHOBAR that will be returned. 
The normally unchanged VBAR 1s adjusted to agree with 
these abbreviated results. 

For NX =+1, supersonic velocity correction 1s operating 
and a message to that effect will appear each time a region 
of local supersonic flow 1s encountered on the body, and 
also when it ends. The body point number where these 
transitions occur will also be printed. 
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ON-BODY POINTS 


Body 1, 2, or 3 


For both compressible and incompressible versions: 


— 


ma K XK 


VCOM 
VBAR 
M ACH 


CP 


RB/RT 


PS/PT 


VINC 
M ACH 


CP 


PS/PT 


The index number of the point 
Axial distance 
Height 


Surface distance from XRI 


Compressible Version 
Resultant velocity with compressibility correction applied 
Average incompressible velocity at a given axial location 


Mach number 





Pi ~ Po 
Compressible pressure coefficient 
Io 
Po/ Pt 
3.5 
: vy 
Static to total pressure ratio, |1 - 0.2{— 
a 
t 


Incompressible Version 


Resultant incompressible velocity 


Mach number (based on incompressible velocity and compres- 
sible flow equations) 


2 
Vi 
Incompressible pressure coefficient |1 - | — 
6 
1/2 pv 
Static to total pressure ratio, 1.0 - —/“ PY 


Pt 


00 
OFF-BODY POINTS 


For both incompressible and compressible versions: 


Rake number 


I 


X 
V 


THETA 


M ACH 
WFRACT 


VX 
VY 
VRE 
VBL 


RB/RT 
PS/PT 


VY 


VRE 


PS/PT 


Number of the point (points without numbers are interpolated 
values at a body) 


Axial location 


Vertical location 


Flow angle, tan7} (VX 
Vx 


Mach number 


Local cumulative weight flow at a given point on a rake divided by 
total weight flow at the rake 


Compressible Version 
Axial velocity corrected for compressibility 


Vertical velocity corrected for compressibility 


Resultant velocity = vx + vy" 


Average incompressible velocity at given axial location 
p/p, 
Same as ON-BODY POINTS 

Incompressible Version 


Incompressible axial velocity 


Incompressible vertical velocity 


Resultant velocity =/\V/ Vx? ss VY" 


Same as ON-BODY POINTS 


i 
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RAKE WEIGHT FLOW DATA 


For each rake the following data are given: 


I Number of rake 
xX Axial location of rake 
IRAK Parameter to describe location of rake on the inlet 


= 1; upstream of nozzle 

= 3; upstream of body 2 between bodies 1 and 3 

= 4; downstream of body 2 highlight between bodies 1 and 2 
= 5; downstream of body 2 highlight between bodies 2 and 3 


WDOT Integrated weight flow for rake 
WDOTCA Specific corrected weight flow at each rake 
MACH One-dimensional Mach number at each rake based on WDOTCA 


Graphic Output 


Standard CALCOMP plots of PS/PT (fig. 9) and Mach number (fig. 10) dis- 
tribution against S/CUTOF1, S/CUTOF2, or S/CUTOFS3 are made for any body 
whose value of cutoff is not zero. 

If VPERIN not equal to zero, then a plot will be made of the velocity flow 
field at the off-body points (see fig. 11). At each off-body point a vector will 
be drawn showing flow angularity and magnitude. Magnitude will be shown by 
relative size of the vectors (a unit vector will have a magnitude of VPERIN). 


DESCRIPTION OF SUBROUTINES 
Figure 16 illustrates the calling relations between the main program and their 
subroutines. 


Program SCIRCL 


(A) MAIN SCIRCL Read all input, call required subroutines for each scsinent 
as requested, plot each segment after points are gener- 
ated by subroutine, list points; test for reworking of 
geometry if required 


(B) STRAIT 
(C) FNSTRH 
(D) FRSTSH 


(E) TEST 


(F) PRELPS 


(G) SUPERC 
(H) FONIS® 


(I) CUBIC 


(J) SIMQ 


(K) LEM 
(L) MIRROR 


(M) XYC ALC 
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Straight Lines 


Generate points on a general straight line segment 
Generate points on final straight segment of a body 


Generate points on first straight segment of a body 


Bisuper Ellipses 
Test superellipse input to see 1f mirroring about y-axis 
1s required 


Mirror superellipse input data about y-axis so that slope 
(1,2) 1s greater than the slope (1, 4) (fig. 12 (a)) 


Generated points on a general bisuperellipse (table I) 


Iterate on input conditions to find bisuperellipse exponents 


Other Curves 


Fit a cubic polynomial between two nonvertical parallel 


lines 


Simultaneous solution of equations to obtain coefficients 
of the cubic polynomial 


Generates points on a general Lemniscate 


Mirror the hub points to obtain the shroud 


Direct Interpolation 


Executive routine for the following modules purpose 1s 
to generate points ''correctly'' spaced along the curve 
defined by a list of input points. Inputs are used to de- 
velop double 3-point interpolating polynomial in suc- 
cessive regions along the curve, Polynomials are then 
used to suggest points, derivatives, etc. which can be 
tested for spacing as defined by standard criterion 
(see comments in SPGEN listing) 


(N) SGEN 
(0) DSTRP 
(P) SPGEN 
Q) DNTRPC 
(R) FNTRP 
(S) FSTRP 
(T) FNTRPC 
(U) FNTRPA 
(V) TLU 

(W) LIMIT 


(X) WPUNCH 
(Y) WRTXY 
(Z) AREAA 


(Z1) DRAW 
(Z2) PLOXIS 


(Z3) SINTP 
(Z4) SORTXY 


ERTRAN 


O3 


Refer to listing for comment card description 


Special Calculations; Output to 24Y 


Generate rake points at requested positions; plot rakes 
Write all X, Y coordinates 


Compute area 


Picture Plotting 


Plot X-Y meridional plane picture of each inlet segment 


Plot frames for inlet picture and label axis 


Utility 
Lagrange three-point interpolation 


Rearrange the values in an array, x, to increase with 
increasing index (ascending order); sort y accordingly 


System Library 


Routine which gives FORTRAN access to several 
UNIV AC 1110 operating system commands. Can be 
eliminated by defining unit 25 before execution 


SYMBOL 
SC ALE 
LINE 
NUMBER 
PIOT 


A) MAIN 24Y 


B) ABFORM 


C) ASSEMB 


D) COMBO 


E) ELFORM 


F) FILES 


G) FLOWS 


H) GEOMCF 


I) MAFORM 


J) MAIN1 


K) MAIN3 
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Standard CALCOMP routines needed for all plotting 


Program 24Y 


Executive calls to MAIN1, SOLVE, and MAIN3 also, this pro- 
gram rewulnds all disc files used, 


This subroutine calculates induced velocities array elements 


in arrays A and B, 


This subroutine reads in normal and tangential onset flows in 
row order and assembles them in column order. Also it 
reads in non-uniform onset flows 1f any, 


This subroutine calculates combination coefficients 1f required. 


This subroutine reads in body coordinates and associate data 
and calculates element data, 


This subroutine assigns temporary file numbers. 


This subrouting calculates individual and combined flow for on- 
body points. It also writes the individual flows on a tape file 
for input to NOZZLEC, 


This subroutine calculates required geometry combination 
coefficients associated with parabolic distribution, 


This subroutine forms and stores matrices A and B and ecalcu- 
lates alpha and circulatory onset velocities. 


This subroutine calls subroutine ELFORM and MAFORM and 


reads in the flow control parameters, 


This subroutine calls COMBO (if required), FLOWS, and 
OFFBOD, 


L) MIS1 


M) OF FBOD 
N) OF FPTS 


O) PRINTG 
P) PRNTEL 
Q) QUASI 


R) RMAX 
S) SOLVE 
T) SOLVIT 


U) TYPE 


V) VXYOFF 


W) VPROFF 


X) WEIGHT 


Y) XYFORM 


(A) MAIN NOZZLEC 
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This subroutine inverts a matrix and Solves simultaneous equa- 
tions. This is used when ISOL = 2 and number of points less 
than 102. If more than 101 points are input program shifts 
to SOLVIT, 


This subroutine calls OFFPTS, VXYOFF, and VPROFF., 


This subroutine reads in the off-body points and prints them 
back out. 


This subroutine writes out the body coordinate data. 
This subroutine writes out the body summary data. 


This subroutine performs direct matrix solution. Program 
uses this to solve matrix when ISOL = 1. 


This subroutine calculates maximum radius or chord, 
This subroutine obtains the sigma solutions from QUASI. 


This subroutine performs direct matrix solution (used when a 
very large number of points have been input). 


This subroutine checks input card type to verify proper data 
card is being input. 
This subroutine calculates the individual and combined solutions 


for each off-body point. 


This subroutine prints out the individual and combined solutions 
(in blocks of 100) for each off-body point. 


This subroutine calculates the vorticity weighting functions for 
each element, 


This subroutine calculates the indiced velocity array elements 
A and B, 


Program NOZZLEC 


Executive calls to INPTR, SEARCH, ANGLEF, SOLVE, 


OFBDY, and if compressibility correction desired 
COMCOR. Also, this program calls PLTER if output 
plots are desired 


(B) INPTR 


(C) CONST 


(D) SEARCH 
(E) ANGLEF 
(F) SURF 
(G) SOLVE 


(H) COMCOR 
(I) ONBODY 


(J) OF BDY 


(J) VBARIT 


(L) INTER 
(M) SORTXY 
(N) INTER2 
(0) INTER3 
(P) INTEG 
QQ) CALTIT 


(R) PLTER 
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Reads input parameters from unit 5. Also reads 24Y 
output coordinates and basic flow solution velocities. 


Calculated most constants and intermediate parameters 
and prints results, 


Finds the highlight on each body. 
Finds the body surface angle for each point on the body. 
Calculates surface distance along body as a function of X. 


Computes linear combination coefficients A, B, C, and D 
to satisfy input flow conditions. Also calculates the in- 
compressible velocities and average velocity V, (used 
in compressibility correction) for each point. 


Applies compressibility correction. 


Uses velocities and densities to calculate: pressure ra- 
tios, Mach number, flow angles and list all results for 
points on the body. 


Same as ONBODY for rakes. Also calculates local frac-~ 
tional weight flows for each rake point. 


Calculates average density ratio for compressibility cor- 
rection. 


See SCIRCL routine SINTP (duplicate) 
See SCIRCL routine SORTXY (uplicate) 
Calls INTER 

Calls INTER 

Performs trapezoidal integration 


Titles plots of pressure and Mach number versus sur- 
face distance, 


Plots pressure and Mach number versus surface distance. 


ov 
INPUT/OUTPUT UNITS 


SCIRCL 2 - Standard card input 
6 - Standard output list 


4 - Temporary storage; if flag J or E is 1, input to direct 
interpolation routines is written here 


17 - Saved. If flag Ais 0, input for 24Y consisting of X, 
Y points is written here 


24Y 0 - Input ©7 from SCIRCL) 
6 - Standard output list 


7 - Saved. Input for NOZZLEC 1s written here (6E13,.8) 
X, Y, V1, V2, V3, ete. 


2-4 
Temporary storage 
8-18 et 
NO ZZLEC 5 ~ Standard card input, flow conditions, etc. 


6 - Standard output list 

7 - Input from 24Y 
12 - Saved. Input data (for body 1) for boundary layer program 
13 - Saved. Input data (for body 3) for boundary layer program 
14 - Temporary storage 


See eee - -_—— - -— ee _ ——~ ~~ ony Se ee ee = _ 
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6 91 C READ IN RHS IN SIG ARRAY SOL VOSC _ 
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@  _=235 #&3=§ Ns ens ¢ NN - eS eee ee ee ee, __SL¥T23©_ __ Ses uaa |. 4h a de 
236 WT = NT © NN SLYT23¢ 
———— 2457 __200 CALL. GETTOMTs 20 Mis AINSDs Le AAZD ees a is, ete ek 
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] 250 W222 ¢ Ww - (1 SLVT25C an 
252 C SLVT252 
@___253____¢ = = wane THE MODIFIED ROW OM TAPE OR CONDENSE THE ROW ———__ fin SUN QS 96 8a Bt Baht rasa tee ee a a 
258 SLYT258& 
258 Ses AU SRV P Ge eee ee ay a a eee 
@ 256 IF ue “GE. NROWD €0 TO 230 SLVT25€ = 
258 NN] = NF CONS OT SLVT258 
@ j|259 ~ 2 US Te ed ES Gi i a ee es 
260 CALL SAVECNOUT, %. NNo NND, AENSD, NN2, AUNLDD SLVT26C 
_ 261 ; 60 TO 250 _ SLYT261 = ea \ueeeeee ee, 
@ 262 230 NF = NL - KOLD SLV1262 é 
268 AINF) = AUCMN) StL ¥T264 
@ i 25 290 NF Pose UNE Oe Nee octyl eet en ees eae a Soe ee SEO 8 yt oe eee eee 2 i 
266 250 CONTSNUE SLV¥T26€ 
—26T REWIND MT oe eee es et. Jt 12 
@ 268 REWIND NOUT SLYT268 J 
270 C + = SWITCH THE TAPES SLY¥T27C 
Ona eo es ree ee es Bessie, tee ee LG OUWN REN rs, coh, ie oe ee ets ces 
272 wT = AT SLYT272 
a 2 OUT _Skvv273. ae ae eee ee 
@ 274 NOUT = NT SLYT274 hee! 
275 C SLVI275 
276 C - = LOOP BACK THRU THE SOLUTION SLVIZ7E 
an Te ges Oe Ree ene ee GE tes Se PT ee ek pce etn ee bl tt 
278 Ni = WF SLVT27E 
pe 2 TO BOL FO 220 ieee Bisie aeM pie sede en ts oe SEVER TS CS ne. es te es ees oF as 
@ 280 C SLYTZAC ~ 
2 - - A UI a es a 
282 c SLY¥T2&2 . 
@ 28326 REWIND NIN ath 8: Son ps. Te eehe  geee eee Le SET 2830 eae 5 oe: eee! ete tJ 
2784 N2 = WN StvTzea ee 
285 C SLVT285 
@ 286 Ct NoTE ne AT IMIS POINT ALL LOCATIONS Atl) _TPRU ALMORE? ARE FREE _ S$LVT2e8¢ _O 
287 Cc es 2A aa ke Font eens eed Satna Ete oe 
































- ? Qo 
291 NS = M1 SLVT292 : 
a 292 i, aod WSN, ete _ Site ust ds Saley prabunlare atte nk SEWI EOS. _ a 
293 Cc St¥1292 
_ 298. = = READIN THE SOLUTIONS = ss _” ot LL SRT 208 ae ims Yeas 
“a 295 C SLVT29° 
a ad rd as a I a Si a Un an 
297 NM =I NB - WS ¢ 1 SLVT297 = 
Mm —_zgR CALA _GETTENINs Lette AINSDs Le AA2D SET ZO ae ; a ne 
299 NT CNT ON SLVT295 
et 2 = LY tre ter certs ete tat Gk fea i sige be a wig al et Sea ae eae oe _SLV¥T30C_. -.... |. _ - a eee, eee 
301 260 N2 = Wi - 1 SLVT301 
303 C - - WRITE THE SOLUTIONS ON TAPE SLV1T302 
® i  j|(.(3s393  ¢_. i ee SLT 5O8 ne pte Page i tee ee A etn 
305 NT = 0 SLVT30£ 
_ NG tC“(i‘“‘C“( ‘( ‘COCNO 290 YO SE 1,6 8 §=—h—h—h) UU i a ce 
a 307 NS CNT eo 2 Stv1T30? 
om C 
309 290 CALL SAVEUNW, Le We No AUINSD > Le AMZ SLVT30S a 
a ABN a a ee we i SI Rn ee we te 
321 CALL TIMEVGAA2D SLY¥T311 
2 BBS AZ = AAD 7 BOW STG pide yo ne ae 
@ 313 WRITE (6.3003 Ny MN, M, 88 StvT3223 ie 
318 300 FORMAT tSHOTHE IS, 2H X 1S, 12H MATRIX WITH IS, 35% RIGHT SIDES MASLVV3I8 
315 1S SOLVED OIRECTLY IN F8.3, 9H MINUTES. ? SLVT325 
@ 6 RETURN Pe Sk a re ah nee ii NO tt rh tee ee ee 
317 ENO Stv7¥3i? 
@ Se er ee ee a or ae eh es Se ee ee po ke ger rs a —_—_—_——— «ir, 
wi 
eee —_ Gas <a ——.— eee eee oe ee as 
— a : 
O______1. —————-SWBROUTINE TIREV AT.) Se cae they oan A RNa, Me eae ER Senne a oe eR | Stray a, RR ees, pte 
T=0.0 
eee ee a ce VAD eee cee — Sa 
© ® Eno 
OW ee Ces oe ee, FS se seeaee es = Sse see eee tt ae Soe etn ee ee 
.— eit dn Meats, «Wiig 5 Rig MOM a ode ae, ano parts 
C ij sic CS RROUTINE YSFIW nee eee i i a 
2 RETURN 7 
Set ca a a Re ar ee ee eee 
EE ooo 
we ene aes fF fess vee tee =. —— 
° a 6 
a. eucle 
1 ________-SMBROUTINE __TvPEtT6000. IBAD? —____ en © 1 J 211115 | _ ok a aoe ee wise 
TYPEON2 
eas ee Sires ele pte hn > tet ete sets PE eee ea Sarge tee onan te ates 
| ww a WRITE (6020) 16000, 18A0 TYPECOS “ 
D HR n ) 2 F P e PF on 
a © CARD. HOWEVER A TYPE *.12, ° CARD WAS FOUND.°/1N0. TYPEOOE 
OF ——2__________<CHECK_ OVER_THE_ INPUT DATA_CARD SEQUFACE FOR °, AYPEGO? ss ss ee, ee eee we Re 
8 °COMPATIRILITY WITH DESIRED OPTIONS.® 3 TYPECOEe 
sisi Sg ed te ae Se, fa & ae ~ oo _OTY PEEOONS ote 7 ee Pe, | eet 
wu 10 Cc TYPEGIC ~ 
a aS RE Og Sg i ED eet oe A 
2 12 30 FORMAT C1H0. “BECAUSE OF THE ABOVE ERROR, THIS RUN IS TERMINATED *)TYPEDI2Z 
- Bree ee a oe AN ats Ss eee. TYPEQL? il 
14 STOP TYPECI8 
Be te ee GR ha Ye te OI ee eee, esas ee ee hee __ ____aYPEQIE eee = ae 





@____1__suprourrne _vavor gen. wo e Xs) OND “ 
¥WxXY¥GOO? 


otro erin, aso, som ESOC 
oVNGRZ3» WVEL2ZI. WNTCIZ)» VYNGI2? wxvooos ” 


6 COMMON/JELD/ NO¢CSOO?, VYOISO0). OS¢S00), SAdSO0d. CAtSOO), ¥Wxvoove 


@. = =phry) cue esnod, pteson gon 
8 ¢c wxvooos 
Ff SAB OncCOMRONsCCh IMCLTCLISALPHASSUMOSILORATLUCIOSIZISIND owavoons 
@ 10 1 o ALPHAO, CNUI10). SMOSWFC10), M1000} WXYOOIC “a 


i2 COMMON /ECF/ COCSCC). CFISOD). CEISNOF, CI¢(SOO), WFtSOO? VN VCOL. 











Dn 8 he 2 ” 
von 
18 COMMON/BFLAG/ TOBt10}, INLCLOD. IFLEIOD, WLE1OD, RIFTCIO)» vrvoois 
eS c = 
® 16 2 BTITLECIO, 7)» I8T, IBST, LBTOTs WELTON, vxvoote - 
i ~eeadda al ry ‘ aga 
e 18 8 “1bha (0). TFST¢10} ISECt1OD~ FTITLEG2S). IPVREIO)  vxvoolE a 
a ~~ Common sconerriit eee 
ae vxYOOIS 


staat neers Ba ae sl ee eee ek 
{ @ vxYO023 “ 
4 24 IFO6o IFN7%. IFORs IFO9e FID, vxyvoo2! 


























Oe fara ret Fs deaee des WO en So Ss ey oe I 
oIFi6. IFI7. IF18. IF 19¢ IF20 vxyvoo2s ts 
oe REWIND IFO} sk Se OOP net et = ees ES 
@ 28 REWIND IFO2 WEYOO27 es) 
29 C vxvooze se ty oe ess 
30 C vxvoo2s 
i 2 eae 31___€ SET sone_ouants Shc. VXY¥OO3C a ee __ 
= as 
re a = - Bes: Ba Sh a) a te OOS ee 
; @ 38 OSI = 0.0 vxvoo32 wet 
35 c xYOO3A ; 
a 7 DO 230 I = 1.NO WXYOOSS 
et Fae eae tnt, ast WX YOO3€ wt 
38 N= 2 = vxvoos7 : 
39 c wxvoos3e —_ 
j oe }§}©=©6—6 5&0 80 TB = TCS TOT << NOOSE — OY 
: ai IF C18MFC(IB) LT. Cd 6O TO 150 ¥XYOOSC . 
a2 C vxyoos) 
@_ 83 ~~ ~CCOUNTER FOR ELEMENT GEOMETRY = Ae ast ai WN N a I Nee ate et ca Se 
ae J = INLCIB) - 21 vryoos3 
as C COUNTERS FOR A.B ARRAYS _ VxXYOOsa estas oe oe eee 
@ a6 JI= Ju ¢ i wrvoces i 
a? JF = JJ © MECIBD ¥XVOOSE 
a8 C vxvyoos? 
{ @ a8 C det TS COUNTER FOR THE CURRENT ELEMENT a OO a a a et ee 
c vxvOORS 
Ojeda ies out ons Bi i we -peebe OUT ACB ARRAYS es - set vonsc a eo fone te oe Rn ees by at ed 
i 00 30 JJl = JTeuF vxvoos! 
tJ = 0.0 oops —( isi‘ SM 
: 58 30 BtyJil = 0.90 wxvyocs: x 
Os SS eet oo Se evens oo 8 tg 
56 JJ = Jl ¥xvooss . 
i@é = " -) 








BR SN I RFI 














57 JJ1l = uJ ¢ 1 VxYCOS6 
ss JJ3_ = Jl + 3} ee Arts ae __WevOOS] ees __ = i 
60 TO SO vxvoose 
7 sae ADN dt IGS GOON es aa ke Oak ea Sete teen ie eae SVE ce eee ee _ a 
61 SaovJ=5s+1 VXYOO6C a 
7 63 1s DLEDDy SALUD» CANUD» ISECCTE)s Aly ACe APs Oe Ble BCe BP) vxvo0s? 
a AaB a ee 
65 AtJJ) = Atus} © AO WXVOO6E 
bh Cc ea! It _—__—_trvonae 
67 Cc vxvo062 
——-___ 4A. dD SIA EO. £) £0 TO ago ln 
69 If «GIF STUIBD - 236C6,70.60 VxXVOU6S = 
71 C FERST TERMS OF FIRST CROER WXYOOT) 
_ 2 tC A OH A DD | Scene ete Ne 
73 Biss ») = Btu Dd - (CFCUP © COCUDDSBL WxYOOT3 
Zs 7S BtyusJl) = BIJs1dD « COiv)e81 wxyoo7Ss 2 
on @ te es Bet eee 
77 BiJu3? = Bisus) ° cFiudep) vxyvoo?? 

i a ah et a Ne =e coke weet Pee ee 
= 79 70 IF CIFSVCIB) sk Te 22 GO TO 8D vVxvOOTS ‘y 
2 81 C CURVATURE TERM OF FIRST ORDER vxvoo81 ss 

_ Addl D = AGS Do OitpecuRnviseac 00 000 dM 
83 Btus } = Biddy Dd © OL CUPOCURVIJIOBC vryoge} i 

a Na oh a 
co as 80 IF ¢ISECtTB) .£€0. Cd 6O TO 110 vxvoos! 

BT ARSE CATR = 239 60 A009 90) 
i 67 vxvoos? 7 
= _ 88 FIRST TERMS OF SECOND ORDER VO GG ai a ee 

69 90 ACJ ) = ASS D - (CEES) © CI¢U) DOAP v¥xvYOO8S 
90 OE BP - CGE Clip aeBP 6 tC“‘“ (“‘ONCOCOCOC*#*#*#*#C*C*C*C*#*C#dNNR VQ Seas 2 
C 91 AtJJI) = AIJUID © C6CUD SAP vxvo09) 
93 AtJJ3) = AtJS3D © CLEUPeaP vxvoo9? 
eS 94 BtJI3) = BiIII3) © CItUIEBP yxyogos ee 
95 Cc VxXYOO9S 
_ 96 100 JF CISECtIB) oLT. 2? GOTO 3110 __ YXYOOSE eee a ae 28, ae 
Oo 97 Cc vxYo09? 
99 ACJS PD = ALS 9D © Ze OCAPSEOL CUI@CURV EU) D002 Vx YOO9S 
~ 00 BY DF ES Bt DD * 2.09BPOTOLCueCURVEUREOOe2 ea ee 
101 Cc WxvO10) 
ig 8 an i i ND a oe es in ae et tl hy ee IO St ee es eee 
103 JJ = JJ * 2 vxvyolo? 
108 CT a Ee 2d 
105 c WXYOIOS 
~ 306 2 IS SE IR MH 2 OE ee 
107 60 TO 50 vxva10) 
46 —__A3Q abel Saba md aaron = 
109 WXVO1OS 
111 IF @LIFT¢IB) .EQ. CP 6O TO 2350 VxXYO12) : 
ge i NO Bas tlk Banke et eh at es I pe ow ee to i ee fete ea 
133 vt = 0.0 vxvoll2 





y — 


+ ‘Mer. 
FEO ey anshetndlshrnansfyae fan n/a chsRaetensnstfy yeu, ersten rere iia = 
1’ 00 1480 J = JILJF WVXV0O114 


1g RK eer ene rer 
114 18D VT = VT © ACJIOWF (UD VXVOILE 
a YS Vn CT ST © 
118 K=K 41 VXYOII6 
; 120 VYNCKD= VN VXYOI2¢ 
a a enters rere ec es aid 
122 c Wx¥O122 


i23 C SET UNIFORM ONSET FLOWS VXVO123 


12s C ALPWA = 0 VxVO124a 2 
125 VXTCJI= Cal VXVOI2E 
“3 126 VYNG19=>-SaTt VNVOIZE 


c ¥0121 a 
128 C ALPHA = 90 wXv0128 


129 WVATE2Z)= SAT eee ZS Fe 
c 130 VYNEI2)= Cal VXYOI3C 
7 132 SET INPUT NON-UNIFORM ONSET FLOUS VXYO132 


ecu eae Ss use he ae he, Sa oe nt ee 
134 IF ¢m FO. Mid 























60 To 170 eyyoiss 
135 M22 MA o } WxYO eek: 
C 136 00 160 K = M2, » pared: 
137 VXT¢M?) = O40 VXVO137 
138 160 VYNEK? = 0.0 VXVOl3¢ 
m 139 170 CONTINUE VEVO13$ 2 SNe fed ents 
180 Cc WXYO1SC 
141 C_ CHECK IF INDIVIOUAL FLOW ws. OESIRED VNVO18) ie se Baek cae sat tenets as 5 oe 
rm 142 If (INO MNEs 1) £60 200 VxXYO182 - 
a3 REWIND IF 18 VxXYO182 i 
148 Cc VxYOIaa = 
a a ASS) CC CCALCULATE INDIVIDUAL FLOWS _____¥xYo1as eats. eS ot Reet cee em. 
Cy 00 190 K = Bem VXYVO18¢é 
187 CALL GETTCIF IG, Le Ne SIGI1D,s Be VND VOT ets ee cee _ ees 
QO 148 WNUKD = VXTOCKD wxyoieg 
149 WYK) = VYNOK) WXYO1SS shee on woe 
150 Cc ¥xX¥0150 
C_ aS? si 280 OM 2A) Se be tas 
152 VXEK) = WXEND © BCCDOeSIGIO D WXVO152 
153 180 V¥(M) = VVEK) © AtueesIgty » __ __ _¥xv0183 aa — 
Cc 15a 190 CONTINUE vxyoisa 
155 Cc VXYOISS a ees 
156 C SAVE VELOCITIES WXYO1SE 
es AST WRITE CEFO2D  CVXERD, VOM, Kt 1.) __¥K¥O0187 nee eee 
158 Cc VxYOIS# 
fee 159 C CALCULATE COMBINED Flow VXYO15S os Sea ee 
C 160 TFUFLLL.EO.0) 60 10 230 ’ 
OO vxc = 6 
162 vYC = oO VxXY0161 
C 163 o9 Te KM = Jom tee VXY¥0162 rasa eer 
168 XC = ¥NC SUNTIKD OK («ed vxvO16: 
16s _ 210 iaceve CS VY¥C CVYNEMDO@CKEMD VXYO164 _ _ Mos iss 
Cc 166 Cc VRYOL6S 
167 DO 220 J = JoN VAYVOR6E 
168 VXC = WNC ¢ BludeCsIGts) ¥VXYO167 6 
169 220 vve_= v¥e + aturecsrctur oeuvre oe hoes 2 ee oe ee 
170 Cc VXYO16S 3 
172 C SAVE VELOCITIES ¥VXYOI7C ~ 
O_o m2 WR ITECTFOLD _¥xce NYG Sed, See i Seem eas VXYOI7TD Be utes Ae a 
173 Cc WxXYO172 
ee ee, | ere Bert es ee ee a a ONO, ae ee ah eet ee 
u 175 230 CONTINUE VXYOI174 e 
€ 
177 c VXYO17E 
~ ——_____ 278 RETURN. ll ts ee ke Lae oe WYOORTID. ~ 


179 Ewo VxYO17@8 ay ae 


Peres El ripe A IO Se emma pct an ence -—— es = > ooo 
a 1 CCMPILERGYM=1) = 
2 SUBRCUTINE VPROFF IN» M, NOg xX. Y. TITLE, IND) VPROCT1 
epee os 3 Cc VPROODZ 7 
a CCMMON /FILEINS IFAL, IFA2, IFO2, ITEFCU, TOS, vPROON? 
Ole TEOk. 1609, Iftn, Rone BiB SS eal es foe: 
6 2 IFC he FF025 YFLZ. IFS, TFS vPROOre _ 
ne Soren Tt og 3 : eIlFl6G— IFL7. IFI@. IFI9, IF20 VPROCOE 
8 DIMENSION XE 9, YC1D. TITLECTD, VEL2UODVINEZOOD.VZNI2000, VPROQONT 
etehi. AiG —~ Lo V3xXC2Z200)eV4XE200 DeVS5Xt2008 : 
10 CCMMCN/OVERS VOXCCCODOVIVEIZ0OdV2YE200 ov tVEZO0 eV 4VI200), 7 
HE a ee pe St ee A ae 
12 COMMON/COP/IFLEL . 
gion, - 12 BS 32. Bis Pe Oea, cee a5 = = Coats : “3 _ wPROGCE 
18 DATA RO/5762957797¢4 vVPROONS 
ess 3 SRS shi, CE Fes Fe a oe, = 5 ak as ey ae _. ___WPROOIC : 
16 C WX AND VY HAVE BEEN SAVED ON UNIT IFO2 ¥vPROO11 
18 REWIND IFO2 VPROO12 
eat Oe i! AO et eat ea te ete ae. og : : nee es VPROOLIS — s 
20 IF C€IND eNEe 13 GO TO 4&0 VPROQIS 
et pa eee Oe a ae i x : _ ; hs a te _ vVPROOIE : 
22 C IWOIVIDUAL FLOWS VPROOI? 
= ¢ PRO a. a 
2a Il =k = M2 VPROO!IS 
— ..@ | 2 S00". VPROO?L 
26 DO 1G Jv = 1.N0 VPROO?)} et 
etait Men Ce ~~. .. Y1 2 71 6 M2_ VPROO?2 S 
28 Iz = 12 ¢ m2 VPROO?3 
9 Cc VPROG24 Ts gate e oe 
30 10 READCIFO2) t¥EKd, KM = T2512) vPROO2¢ 
Seeiet des, al Cc eats Toe = VPROOZE | 
32 DO 30 K = 34M VPROG27 
a 33 ; WRITE(606C) mM, TITLE ~ : _ .  WPROCZE 
34 C VPROG2S 
35 O00 20 3 = 1,NO VWPROO3C 
36 IX = ¢]-))0"%2 + 2eKn - j VPROO3) 
7 37 Ivy = Ix + } ; VPROOtZ 
ai 38 VT = SURTOVEIXD82 © VETY) e082) VPROOX2 
39 TH = ATANZEVOTY), WEIN) D2PH VPROOIA 
40 FFUM EOL) VIKCTDoVEIXD 
4] : JFUm ECOL? VIVETI-VILYD Soe ae 2a Batata cats 88 fa ee 
42 IFIKEGe2) V2NETSZNEIXD 
ee 43 IF CM .EQe2) V2VETIIMEIY) 
yq IF(K.EQ.3) V3XCITD COVED 
oo bees 4S IFOM EC C3) V3VAITIIVEIY?D 
46 TFUK ECLA) VUNETI=VEIXD 
47 FAK 4 eVtIo- Soh a OE Mae? io ae 
9g TF¢CK.EC.S) VSedToprvelxD 
a ao _ IFW .EG65) VSYITIZWEIYD _ 
50 IFUK FO eb) VEXETIIZVEIXD 
eee 51 IFiK.EC.6) VoveTpovely) 7 OF 
52 C VPROOSS 
53 _2O_WPITECGs7C) Ty XCIDs VRID), ViIv), VaTYD, VT, TH __.__ VP ROQTE Bob tek es sg as ota, een 
54 C VPROO3? 2 
ae 55 3¢ CONTINUF VPROCTF oe 


56 WRITECTSIEUO) CXb Ube dZL NOD ° 


Se ee ee ae 


eo 690008ee60388626¢6D80)hlrtC~—~MSDSC(“(‘i‘i‘i‘ DS 


arene by 


“~ 


“_____ 





-_ Pp 


57 URITEC7.1500) (¥¢J)e0=3 NO? 
“______.38 _____write C7olSQO) UVINGDed=teNOD ee, e 


a RET OO) Wx t3P oJ=1 NOD 


a a re 


WRITEC7, 2500) C¥2X CI) eZ NOD ee eh oe 





URITECT 1500) CV9NCI) od=2 NO? a 


63 I 


FaN.GEo6) VRITE(7 250) VEXED JH? NOD 


pee Terre e1SON) AVIV Gddot=DMOd 


VRITEC7T.1500) (VZV 42) -3=1 NO? 


eee ae ——__—-WRETECTsISON? £WSVMdRIHI MOD 


WRITECT,1S00) (V8YtJs od=1 ¥NO) er 
_... ma RSL FORA YAP GE Dg a OS 
71 Cc vVPROO3S 


Fa ____ 2 menen 
C COMBINED FLOW 


| 
| 








73 VPROOS } —" 
75 80 WRITE&6,80) TITLE VPROOS 2 
gh ee a et et ee 
TT 00 $60 1 = 1,640 VPROQSa 

READ UTBOND tad, ¥ Roose 
79 VT = SORTEVELI©e2 «© vV«I2)082) VPROOSE Z 
8) SO WRITECG TOD To XCIDs VEIDe WELD, VEZPy WT, TH vProose aa a 
a2. ROS Oise. see ee ee A. eves 
83 60 FORMAT(IH1, 20K °TMIVIOUAL FLOW NO. *oIZ- SK*OFFRCOY POINTS*, vVPROCSC pas. 
6a L__SXsTAS¢/T5 21" s 127s HET) %a TOs VELD's T6LeVE%, 9 wPROUS) So Mose heta eet. enh es, 4. 
85 TTB *V¥*, TUHe*VI*, FIND, *THE TAIDEGD* 74) ¥VPROOS<. 
87 80 FORMATI1H},. 20x “COMBINED FLOW FCR OFFBODY POINTS*, SX,TASSss vPROOgs: a _ 
68 1 Tas,° I*® P T27,°RCTI*%, F T948,° *vt(rI5*, T6R,*¥X*, VPROOSS See See ew 
89 2 T7O.°VY¥*, TIS oT 11096 THE TACDEGD°7/} vPROOSE 

sc a a a ee he ROOST ele See 
9) 91 RETURN vPROOS® 


22D poneS 


kK 








_ = mn a gn a ge pg are, cat. == — —— a a a a a es Se, rere 























& 
ee oe anmmn nner neem ne nee ee A A I A LY A ce ge 
@ a FUNCTION XYFORM. CXLA_VIe OST SINT, COST, IFSTe = vom i 
2 2 MJe Yoo DSdy SINJe COSS, ISEC, XYFOO0: 
EO ae tee 2 A AC e AP» BO, Ble BCe BPP ss VF OO0 Se ee 8 ee eee 
& & c xYFOO0S 
€t 
6 C THIS ROUTINE ACTUALLY CALCULATES THE INDUCED VELOCITY XVFOGOE 
& 7s ARRAY ELEMENTS _A-AND Be i esti dm sass = 
8 Cc xvFoode 
ye Sa a PA aed at ls hk Te ta a, A alee ee Naa es Se eS OE 6 ee i te 
10 DATA El E2/169eDel lo LLISeEOeEY/0 2000 1600017 XYFOOIL A 
ee beet es > I an ri ee dt 
12 c XYFOOLZ 
eo sa ON St ule ne oo _ OND oo ee ee espe oye et 
14 OY = v1 - YJ KYFOOIA 
_i8 ss  ROSO = ONO? © Pyeer ois Se ee ee et ee ee 
16 DSJSQ = DSvee2 NYFOO16 "= 
18 IF €ROSO eLTe OSUSG#E) GO TO IC KYFODIE 
®@ 19 Sa rath vee ONE Oc ry es Paes: 
20 C USE FAR FIELO FORMULAS KYFOO2C 
2i).0Ci OS Ze OOS IZRPOSQ VF OOZL Be ead ees ee 
S 22 VY = vxeDY XYFO022 *y 
28 AO = -¥NOSINI ¢ VYeCOS! NYFOO28 _ 
® 25 CORN = AO eo a eS | |. ee ee er ree 
26 @0 = VX*COS] ¢ VYSSINI RYFOOZ€E © 
27 _____ EF CALF ST 60£G- OF ANDe SISEC -EQe OFF RETURN  —  xyvFOQ2?. vey 
@ 28 $C = SINIECOSS KYFOC28 
9 = SINI®SIN oes 
30 CS = COSI*SING RYFOCTE 
eo. ok CC = costecoss = MYFOO3) ee 
32 $SCC = €88 - COPOVXOVY/12.0 KYFOOS2 
re te a ee Cssc = 4CS * SCHOVNOVYsI2.0  — ZVFOO3? ees eR Ae ee ieee 
34 OXR = 120 - 2-0#DX082/R0S0 NYFOO38 4 
= - sDye - - 
36 DOR = OSISQ/4IROSO9¢.0) KYVYFOCTE 
& 37 c ied 28. Ses, e ___MYFOO37__| eee Sr 
38 Ai = OORStSCeoxR - CSeDYR? - $sEC xYFOOT%E 
=e 39 B) s-DORSICCeOXR «¢ SSeDYR) * CSSC — vFooss. ——si—ss ey oe ee 
@ a0 AC =-DOR*(SSeDxR + CCeDYR) - CSSC XYFOOSO s 
= Re + = s = ¢ XYFODS } 
a2 AP = £07120 XYFOOS8z 
e a3 BP = BOsN2-0-—————‘“‘—s—‘“—stsS ein ooo ee eee as MYFOC€83) is no ert,  cdewes, ee 
aa RETURN XYFOC@A 
. 8 ae ae hile ee ee i rope LL ace el ae 
@ a6 C USE NEAR FIELD FORMULAS KYFOOSE ~~ 
a7 = 9 ° ys & x¥ yg 
es Y = OYeCOSY ~ OXeSING KYFOCSE 
® a9 S = SINE®COSy - COST#SINS _ YF OgRS ea eee, ee oo LL 
So C = COSTSCOSU * STRI*SING XYFOOSE 
som Dh So ao ah ee ssa eae ee et de i eee ee, ee oe ee OOO ee ea es 7 aS ds _— bro uy Sees 
@ 52 IF ¢ROSO .GT. DSUSC#E2) GO TO 2C XYFOOS: Ww 
2 
Sa C USE EXACT FORBULAS xY¥FOOS5a 6 
@ se . _¥B = xX/0SJ st _ , =. ie AO, a KYFCOSE eo _ eA “\]* 
56 YB = Y/DSJ XYFOCSE . 





ee ke 


ar ale ei oe 











RTE a a  -S  e g  e 


57 RO = XEee2 ¢ YRee> XYF OCS? 
eo sa. __RO_= XBee2 = yRee2 _ Bet eee ens as es, 2 _ -KYFOOSE 
59 R1 = RO * XB ¢ 0.25 KYFOCSS 
= -) | -- R2_=_ RQ - XP_* Me25 _ _ Ce ice eee | ees. eat ee. eee VE OCSE ee Ses 
& 61 C XYFOC6) 
thls we tant A? wc = © 4 * = 4 € - 
63 Y = yveosu XYF O06? 
@ _ oa__ X_=_ROSQ - 0e25*DSuSQ_ : are _ ORY 0068 ce es 2 sce = 
65 VV = 2.M@ATAN2ZEY,X) XYFOO6* 

Sata os __IF_AtIFST_e€Qe O) aANDs CISEC eEQe 09) GC TO 30. °& 4  £_ _xYFaQ6E ae a ieee, eye een ae 
6 67 C XYFOC6? - 
sch 8 ne ee es ee 

69 VYO = VY XYFOO6S 
oe mm = + - OnoTrL Sachs ~ 
71 VY1 = XASVYO - YBeVxO KYFOO?) > 
Tike __eWKC = 2 NOVY) © RBOYBsIRIOR2 €0072 Sees yaa 
@ 73 VYC = 2e0CUVX1 © 1.09-2e8tROSO2 - De 2S@RDIZERIORZ? XYFOO?3 ~ 
ss a = = 
75 VYP = RD@VYO - 2.CeYROIXBeVXO - 1.0) MYFOOT7E 
@ mw  t& Bot Sing ae ees AS tal ROY . 2 3 ee beatae «eee ia, eee 
77 60 TO 30 NYFOO71 
Ss sn StS PON On et ee ce ert hee ie ag pe ee an 
@ 79 c RYFOO7S ~ 
Nese ia Se SS a 
81 20 AE = XeDSJ/ROSO XYFOOB! 
ee 82 =. __ B8€ = YeDsusRose 2 NYFCOSz Uther eee | 
83 ASQ = X*82/R0S0 xVvFOOR? — 
a _E€SQ =DSUSO/ROSO __ a MYFCORQ® pee ke : & 
& 85 VK = 2.0FAF@tL.0 ¢ CASO - 0.675) 9€S0/3-0) KYFOCAE ~ 
86 oe Y = eRfs ¢ a eESOs KYF OCRE a 
87 IF CUIFST £0. O) cANDe CISEC c€ Qe OFF 6O TO 30 NYFOOR? 
@...--: 88 WG = DelS*tM.OFASQStASO — 160) % 1.098ESQO —vroose ——esssia—isi—its mae ie oad 
89 WX1 = E€$00¢2.08ASG - 1.0 * A8D76.0 RYFOORS 
OO, _WV1 = AESBES(1.0 * Oo3#tZ2e0@ASO - LeOPPFSOR/3.0  — — = mvrooee sis ene ; 
e 91 VIC = AE*BES(1.0 * Oo3etbe0@ASO - 16098ES0)/3.0 RYFOC9) ot 
_ 92 ‘. BSO = Yee?sROS XVFOO9 
93 VYC = ES00(2c08BS0 - 160 -0.59A89/6.0 XYFO093 
e@ om VYP = @160 © €ASQ - O-259ES000.69/6.0 _ oa _ YF O09 a le 
95 WxXP = AE SVYP NYFOO9S 
96 VYP = BESVYP sees LF ODF fe ee ae oe geen whee oO ae eee 
@ 97 Cc XYFOO9? ot 
_ 9B xvFOOo8 
99 40 CONTINUE XYFOO9S 
€ 300, XYFORM = -SevVxX ¢ CevyY vise VR ONOE ce Joon ok ey ee ee 
102 BO = C8VX ¢ SevY XYFOIO) 
102 __ TF CUIFST o£ Oe OF eANDe CISEC o£ Ce OF) RETURN  —  _____— xYFOIOZ __ me ee Stet: Gwtae tees Stet 
¢ 103 Al = -Sevx) * CevyYl XYFO102 w 
_ 2£oe ee Oy Oe ag i a a se eee 
105 AP = -SeVxXP © CevYP XYFOLOS 
C 106 teers Me, Se Seem bee sheets a: spent 28 _MYFOIQ¢ So eS. 
10) Bl C#VXl * SeVYl KYFO10? 
10e BC = Ce¥xC + SevYC Se — enroicnt a tes Se aay ot 
¢ 109 BP = (OVXP © SevYP XYFO109 we 
110 ( * _ MVE OIG) one 2k ae he On et ts 
111 NYFOI1! é 
a 112 RETURN es, ONDE _ _ _ oS 0 hjs sec es __cte]> 
13 END XYFO112 : 


— 


ee bea 


- 


ly ata = 


ee ee eee) 


—_——— 





SIN SANG peep Ces = a2 ee es 
Program NOZZLEC c 





















@ 
ae | COMMON /COORDS XON IT ON XOF F YOFFt ’ ) = 
2 1 $14700),XTEST NTE STI gXTEST2Z,VCL -YCU eYCL 2, VCL2,YCUL, 
3 VEU eg MRD eo XRZ eK RH ye VRE e VRZ e YAH CUTOFF! CUTOFZCUTOFH, FF Fs 
& 4 3 ELNO,ANG E700) ,ARC 700) -AROFF (200) 7 
ms O 5 afe nn 3 AO if , ro 
6 1 V2"(Z200) pV 3XtZ00) VOX 1200) »¥SXEZOO? pVEYEZ002, 
& 7 2 V2V 1200) »¥V3¥¢ 200) ,VaY (200) »¥SV¥1200) a 
8 COMMON SCONT/ VC eVE1 eVS2eXMC oXMC 1,XMC2,WOOTC,UDOTC 1,b00TC2. 
9 1 TITLE (3) VINE DALF & Ao B eC aD pALC 0A AIC gASC pASCSA 
@ 10 2 B21 ASL oAV1L AST AL2 A229 AZZ eANZ gp AS2 WIC oVICI,VIC2 ~ 
11 COMMON /COUT/ NT ONS! NH NP oI ONX gKND »ICOMP oK ONNHIL NAHIZ oNXHIS 
12 1 NSTZeNSTINST7 NPFRESOD,TRAKC3SCD M1 MZeICOMPI IPL « 
ee 2 HUB a 
18 COMMON /CONDIT/ TICTAL sPT ePSTAT,ISTAT PSTATCeATOTAL ePIC,RHOST,. 
NB 1 RHO TOT sASTAT »QCINE pRSORTC 
& 16 COMMON WRITS AALCARZCeAASC gARSC AASC eAAIL SAAZ TL eARIT SAASI,AASI. = 
17 1 AA12,AAZ2, AAS 2 pAANZ PAASZ 
18 COMMON /SOLUT/ VBARE TOO), VBAROUZCOD VINCE 700) VXINC E200) » 
@ Ye 1 VY IAC (200) RHOBIPCCD»RBORT 4700) RHOEC E2003 7 
20 2 VCOMNt700) .RBOOTI2CO), VRE C2CO), VRECOP E2003, 
ee 2) 3 VXCCMH 4200) ,¥VYCOM (200) , THETA(200),PSCPICC7O0) CAA chee 
8 22 4 PSOPT E700) ,CMACHE 700) .XMACHI700),CP 10700) CPC (700), ™~ 
23 5 RHO1¢700)9 
28 COMMON SSOLUTO/S PSOFPC(200),PSOFPI20C),CHACOC 2009, HPACOU 200), 
@ 25 1 R¥CO1¢2003 ie 
ee ee OR eg NOXX XMIN, —— 
27 COMMON /CLPLOT/ XPEN SYPEN »NXG OMY SIPEN LABEL CIO) -YLABELE 103 = 
oe # }2@©»©»©=6—8 ~ COMHG KT p JCHK 2 LICHK eo eee eras - . ee 
29 C ecween nn ceeea------ www ce a oa eo na aw we nw eee ene rene sees ee cn se nase - 
30 Cc 
@ 33 Cc ~ 
3 C Wars YS we WAYN PROGKAM VHICW CALLS THE SUBROUTINES To 
33 C COMPUTE THE 2-D COMBINATION SOLUTIONS FOR BOTH COMPRESSIBLE 
e- yu” CC ND INTORPPESSTOLE VERSIONS. 9 
35 Cc 
36 Cc 
@ 37 DIMENSION XPLOTETCCE pYPLOTI7O0) »WKKC89,P C19) J 
SG OM TMRTR CS re oP Se ne Sw my eae eT eee ee ee eae ee a 
39 AAAA=VINF 
@ |: = A 
4) IFCICOMPL.EQO.1) ICCHP=0 
€2 TFCICOMP1.EQ.0} ICCHP=1 
@ 43 IFCICOMP ZEQe1) VIAFAVINF OED C—O. 240 VINFSATOTAL F992 982.5 J 
— 7 it $ eet Sy cee ee mF Say ee te ey ee oF ee 
sl CALL ANGLEF 
oe” - =COC—ti‘“‘“SNN ~ WESRVESVE er J 
at VCAS=VS1 mira nna 
as ¥C2S-VS2 
eo 9 a ahs WC WU et, pee ee ee w 
so VSi=VIC! 
che, gon, ets _ _WS2=VIC2 
© 52 10 CALE SOLVE a ee ee ge ot Seeoa 2. « wd 
Me hye, | ee «WC = VCSAVE 
58 VS1=VC15 o ‘ 
eo. 55 VSZ2=VCZ2S = ie See eee J “de 
56 VINE P~ VINE — = — . gee eg ae, oe a > . ’ 
re - aS GS SO, ee, ae eee ae eek ey iy eT ee —_ “ 


57 VINF=VSAVE . 


oN 58 IFtIHUB.EO.0) GO TC 13 : 


$9 WRITE €6,1000 
60 IFUMLOEQo1) WRITE (E2105) ALD AZT AST ADT OAS20A129AZ2,A32¢ _ 
= 61 1 AS2,A52 









6 LF Gm WRIT D2 
63 1 AG2,852 
2 68 ——_LFLUM1 E000) ce AND 172060009) WRITE (60115) AJC QA2C AIC ARC ASC = 
1 


65 ALL eA2l A312 ,091,051 
66 WRITE 1120) We BoC WD eVINEP 

cs 67 GO To 12 oS 
68 11 WRITE (6,10)) 

“a 69 WRITE CGe1029 ALCeA2C eASCARC 


7 WRITEC ) VINEFP ai 


71 12 IFCICOMP.EQ.1) CALL COMCOR 

72 _CALL ONBODYy : eee 
os 73 CALL OFROY se 

74% Cc $$ IND NACH MINADR WO eee 
A 75 C THE FOLLOWING CODE PLOTS PS/PT WS $ AND MACH NUMBER VE Se 


76 Cc on 

77 IF(CUTOF1.LE.0.0) 60 TO 20 

78 DO 15 I=3 NS] Se 
oo 79 IFESUTD RTO) S2ETD=-SUTSSCUTOFI os 

80 id = 

61 15 CONTINUE oul [0A EE emmnmenreecreeeeee 


82 20 = 
Ta oo eo, J1=0 


O00 26 I=1,N$1 





co | oe 
st) 








Cc IFCSAtCI).6T.1-0) 66 TO 16 

86 IFUL LT JCHKID Uryel a a 

87 IF ¢JeEQe2} Sz] 7 
CO 88 IFC] GE eJCHHID JZ 10} Re es 

BT CONVENE eee 

90 IPL=} See 
© 91 CALL PLTERTI1.JeICM1elSei) 

92 20 IFtCUTOFH.LE 00.0) 60 T0 30 : es 

93 IN=NS] 02 
© 94 1T=NH 

95 DO 25 I=ta.¥i 


96 TFUS¢CTD LT.0.0) S1ET3=-SUED/CUTOFEH 
97 1FCSUI).GE.0.0) Sitii= SCTI7CUTOFH 
98 25 CONTINUE 
99 J=0 
O 100 J1=0 3 
“106 0 26 I=Inw, 
102 IF¢S2C13.GT21.20)3 6C TO 26 . 
YQ roy" a (TLE .JCHK 29 Jive] 


108 IFIVeEQe1) IS=] ra 
105 IFCT -6GTeUCHH 2) Jisul ol 
0) 106 26 CONTINUE epeccmeeg 
107 TPL = = CT —— 
108 CALL PLTERIUI SCHR 2,158,123 Sa Ned oe ee 
=e @ oa 109 30 IF ICUTOF2.LE.0.0) €0 TO 37 
110 IN=NH¢ 1 
111 O00 35 I=IN,NT é 


= Ssitn, Ade? TEUS¢T dL TOO) SIITI=Z-SETISCUTOF?2 oS 
113 IFES¢Y).GE.0.0) SLtUTd= SEtTI/CUTOF? : 
sey a eee ee hin ca oe ne a he 


—_ 
v 

















Ogiten Poy eee Ese s . We eS ie Se eR ee iN a eae Ry hs ls in oie eek athe ee 
P Ls sae an ee ee 2 CIYNOXIE TIT 92 01dx a — Ofte 
°° > 0 wi | te TeLLI=1tl 691 e@ 
; TeLI=IL S€2 391 
Se ee eg IHZ*SES*SEZS INIAA-QUOGAA-CIINOAD SE OFZ 491 

o heya 08 es OZ OST *OEZ CNT MK dA IMDOKX-CTINOMDAT == TCG 
So en a _— 62, _ a OhZ OL 09 ENIWATLA IS CIONOAD SI S91 
OeZ Of O09 CNINK* ET CLONOKDOST e9T 

ee ge et PR Sy ee se - HN ®NN=I Ow2 00 £91 & 
LeTSN=NN zt 
OS2 34 09 CO*D 2° GNHIIAT 191 

r OrI3 oat 6 
ee See SS SSS ESS SS SS SE SSS Se ee II=(99 "NN 6SI 
es IANIGNOD 002 gst 

Ba SINGS yee ga gE eee CIYNOAZCIITILOIGA Ls! @ 
CLINOXZOTIIDL O1dX 951 
———— 0 SS Tel1II=-IIk SSI 
c : _— TeIQzI] SEt aSt 
eg ee hE — ee ee ee or OIZ*SEL*SEL CNIWA-OQHOSAA-CIONOAISI OFT EST 

é 


B°0 = (6d 
aie — guOeAL owtwa = (Ltd 6£1 
NIWA = (99d ert 
a ea ee ae Pe } AA = (504 “gt 









” Oh Of 09 (O° O° NINIGANSE 921 
<< $ $$$ $$ ._ Org S TAY IE — Set 
3 CZ2°SI* EMAIL EP SC VasA Td VIO eZ1 we 
aS es nn ””!”~”:CC os 
JNNIANOD OF Z21 
= et =e CoMHIO "19° LDSE 121 we 
I=SE (T9003 dSI O71 
¢Ffz wi } 6T 
a. 9f OL 39 CO*CASD*EEIISIAT SIT 
lll NN YD OF 0O0G00”~™—™”™~™~CAT 
oztr ott 
=f otf w 
7 JNNILNOD SE wt 
Sec Nses Yale Net tet Ske tS ss a ee ee ee eee 
cate oe a Ae aches oe eae — ww 





2 fe ee 





© Bea aie eee ~?} 
te 

a 
171 YPLOTCIII}=YONIID 

@ 172 280 CONTINUE a a pee ee . 
173 MAKKGC7. = TZ 
178 250 11-0 $e 

@ 175 NN=NHe) o 
176 00 270 I=NNeNT eee 
177 IFCXONGID CLT oXMIND 60 TO 270 

& 178 ____IFCVONUTD LT OYMIND 60 40 270 a 
179 IF QXONGTI-XXSEXEP- IND 260,260,270 
180 260 IFCYONCI)-YYCORD-YRIN) 265,265.27 

@ 161 265 II=31¢1 “~ 
182 EI=11}e} —_— $e 
183 XPLOTCEIID=NONCI D 

g 188 YPLOTCIIID=YON(T) “a 
185 270 CONTINUE 
186 IFCINUBCEQ.O) ANN CTD=ATE 

eS 187 IFC IHUB ONE .OD KKK ¢C8)=21 ~ 
ise po 260 I=1,3 
189 NLABELCID=TITLEC1) 


@ 190 260 CONTINUE “a 
191 XPEN-0.0 
192 YPEN=0.0 

@ 193 IPEN=<+3 ~ 
198 NX6=-38 fae aes 
195 NY=0 

& 196 CALL CALPLT¢XPLOT. PLOT KKK,P? _™ 

He ee OO SOP SS 


198 100 FORMATE? 26K, 9 V2 IZK oP VZ oD AKo VI LING VAT LIANG HSS) 


L6I 

























& 199 LOL FORMATO 016K 6 © V1, 12K, V2 oh 2Ky U3 gl2hy VA o/s) | 
200 102 FORMAT(2X,°CONTROL °s 3¥ oSUIPE 10.3 4X) ) 
0 OS FORMATES 10K 6 °A% gd IX, Be TSK» C* ole WINFP 5 / o6Xe80IPEL0Uc3 04ND) 
@_ 202 105 FORMATCE2ZK»*LOWER® 492K, PASSAGE ® p 3X pS C1PE 1003 cA De sy 2p" UPPER" oS, % 
203 Y 2X, *PASShee osx SUiPe 10.3 oaxd dD 
208 11Q FORMAT (2K, °CONTROL %, 3K, S50 1PE 1003 8X) 9/4 2K "UPPER "5 452K 9° PASSAGE" g 
© Xe O.t,9xD) eo 
206 115 FORMAT(2X CONTROL °, 3X p50 1PE 10.3 8X) oS eo 2Xs*LOWER® oJ 02K PASSAGE ® 
1 3Xe5¢(1PE10.3,9X)) 
€ 208 120 FORMATE 10%, °A% oe LIX, BoD IN o°C* 13K e Oo IIK, Se 
—O”””””CSN BX STIPE IO 38d) = 
210 END 
ge OE ES ee ee 
EERSTE RN a ; aL a reper ete 
@ wd 
a 
e- eee oe ae pee ee a eee 
OL ge 
& ws 
SS Ih Fg et gt a ee te 
rn ee Pe ee ee ee a en ae ee 
- ee tn te i NS og ey 
nee eae eT; ee : _ ae 


Ae tet oe 






















en RT NH PN eR Ee 
sae ae es 
om 1 SUBROUTINE INPTR 7 
2 COMMON /COORD/ XON (700) eYONC 700) ,XOFF 4200) ,YOFF(20@3eS¢700), 
ead 3 1 $11 7009 NTEST »XTESTI sXTEST2Z YCL oY FCL1,YCL2,YCU 
a ® 2 YOU Ze XR 1 gMR2— XRHe VP], VRZ op YRHSCUTOFS .CUTOFZ,CUTOFH, = 
5 bk R 
6 COMMON /VELOC/ V1t700),¥2(700) o¥267002,¥40700) oV5¢ 700) pVIXtZ003, 
& 7 1 ¥V2N(200) p¥3xXE200) -V8X4200) -VSXI200) ,VIVt200) A 
& 2 VZY(200) »V3¥( 200? .vV4¥ 1200) -¥5V(2003 
9 COMMON SCONTS VC pVSI o¥S2—XMC NMC LANMC2 WOOTC, UDOT WOOT 
a] 10 i TITLE 03) VINE , ALF AW eB oC oD gALC gA2ZC AIC ARC yASCLAIL 1, “ 
AZLsA ASL eA A32,AS2Z,A 
12 COMMON SCOUTS/ NTN S1 oNH NP oI wt oNX oKND gICOMP OM pNXHIL NXHIEZ NIN] 3, 
Oo 3s 1 NSTZeNSTINST? BNPFRE3S0 De TRAKE ION ML MP 2eICOMP] LIP ~ 
18 2 IHUB 
15 COMMON /CONDITs TI PT oPSTAT SISTA oP T t ee 
& 16 1 RHO TOTASTAT,OCINE oRSORTC 
i7 COMMON/PICT/VPERIN aXXo XMIN GE XEP «VV, VMIN cORDSEMSTOR oAL pAAAA 
18 C wenn ne ow on nnn +e so eee +e On a enn cn en nn ce ew een eee 
* 19 C - 
20 C THIS SUBROUTINE READS DATA FROM BOTE CARDS AND DISH FILES. 
21 c 
@ 22 READ(S.100) TITLE 2 
23 WRITE(6,1013) TITLE 
28 READI5.110) NTeNS lo NHeNP olWoNX MND, ICOMP I ,PHUS 
@ 25 VRITECGe111)9 NT ONS 1 NH ONP ol ONXSKND, ICOMP1,6IHUB ze ieee ee begs 
— £60 ° Ce e¥ ° WP ALF ALXAC, Cie Cle ale = 
2? WRITECG eR 2ID VCoVS eo VSZeVINE cALFAGKMC XMCIOMMCZOTICTAL oPT : 
@ eiéO ebDOTC, ow * s ® ® e ms 
29 1 CUTOFK ,VPERIN 
0 OFT .CUTOF2. 
e 31 1 CUTOFF .VPERIN 7 =F 
—”—«dASC Fe A WPERIN NE 0) READ CTS TSO) MM OKPINGEREP SO VV YMINGOPD 2000 — 
33 IFUVPERIN.NE 00) WRITECE 1319 RX ONMIN,EXEP ,YY,YMIN,CPO 
re eg ——— a 
35 C NWT = TOTAL NUMBER OF CN-BOCY POINTS. oe, 
36 C NSIT = TOTAL NUMBER OF ON-BODY POINTS ON BODY 1 
2 37 C WH = TOTAL NUMBER OF CN-BOCY POINTS ON BOOTIES 1 AND 2 — ’ 
SG CMR COT = OO - «€WEYGHT FLOW CATA AF CONTROL STATIONS INPUT 
39 Cc = 1 - MACH NUMBERS AT CONTROL STATIONS INPUT 
O70 TES = VEL OCL UTES EV CON TROT STATIONS TRB 3 
41 C WX = 1 SUPERSONIC VELCCITY CORRECTION APPLIES. 
i C WKND - DATA NOT SCALED. 
& 93 Cc = 6 DATA SCALED 8Y CONTROL STATION PASSAGE HEIGHT. J 
as — € YCoRPYL = 0 COMPRESSIALT ¥ 
a5 C = 1 INCOMPRESSIBLE VERSION 
oe  }§}3=SCSC C NOTE -- Two H PREE CONTROL STATIONS® DATA mUST EE INPUT, IF Trup N 
47 C IHUB = 0 NO BOOY 2 
ab Cc 
@ a9 c 2 
50 REAOGIS.180) KTEST,YCL.YCU 
51 WRITECG,191) XTEST,YCL,YCU 
e 52 “TRCTWUG EGO) GO TOS 
53 READ(S,190) XTEST1,YCL1,YCUL 
54 WRITECG AGL) XTESTI,YCL1,¥CUl ae 
@ 55 _—_ REAOUS e180) MTESTZe¥CL2eVOU2 ee ———s—‘“‘“‘“‘“‘“‘“‘i‘i‘i<itstsstssss ee fe 
CG WRITE CO LOL? KTESTce ¥CL2,YCU2 . 





ee re ee == oO ee 


6S ees elle eS et 


eat om 


ow ~— ae oa, wee we 
@ 
NS 


® °®@ 
0 
& 


“re 





2 | aA ST A er A Nr ee —w= = = eo ee een 

a Ee ee ee = CO 
5? 3 READ(5.180) XR1,IR2—XRH 

om 58 WRITE COo2 S12? = XR1,xXR2,XRH ~ 
59 READ (S.180) VERE YRH ep ee 


° 
60 WRITE(6,14)}) RJ »YR2,VRH 
a 61 RTES TAXTE STL yxTE STE ARE THE eo OF THE CONTROL STATIONS “ 


c 
62 C_ OOWNSTR F Bo T AND 
C AND STTUEER BODIES 2 AND 3, SESPECT NEES. 


& 7 6a ____ _€ XR1yXRHeXR2 ARE THE PCINTS OW EACH CF THE BODIES WHERE = 
65 C SURFACE DISTANCE EQUALS ZERO. 
66 =, READG7,150) (KONG De J=] NT) 
& 67 READ(7,1SO) CYONEUD,UZ=1L NT) = 
68 READ (752150) CV1Es) w= NT) 
69 READE7T.ISOD CV200) JT, NTD 
& 70 READI7015O) CVSS) pJ=1,NT) “ 
71 READI7,150) ¢¥VRtu) ZT eNT? 
72 READ(7,ISO) CVStI) w=1,NT? eee 
73 IFtNP.6T-100) GO TC 8 — 
7a READC7e250) CXOFF I.) eJ=LeNP? 
75 READ(7.150) CYOFF OL) ,J=1,NP) 
7% READC7,150) @ VIXds) pJ=1,NP) ~ 
77 READG7T.2150) € V2X tu) eJ=]_NPD 
78 READC7ISOD € V3XE.) pJ=2,NP? 
READ(7,150) € VONILP,J=1_NP?D “ 
80 READI7,ISOD ¢ V5X4.) J=JQNP) 
al READI7,150) € VIV¢ CP eu=1,NP) 
READ(7.150) &€ V2V¥ 0.3 ,JS22,NP) a 
OS REDD UF ASOD UC V3V On) ed=leNP) ee me 
ry) READ(7.150) € VOY Cad eJ=1,NPD Sr =, 
oe  —COsS TT ET SO) VCS S=ieNP > T 5 oe 
86 60 To § at 
a 8  @ READtIT.I50) (XOFF Cu) ,J=1,100) 
READE7,150) «YOFF 4.929.321.1003 % 
AD(7, € Vix¢sd.J=-1,.100) 
90 READG7,150) @ W2X yu) eJ21-100) 
OT RERO CF ISO) 0 VIN bu) pJ=10100) ” = a SS be ee 8 7 7 3 
92 READ(7,250) ¢€ V&Xts) oJ=1,100) et 
93 READG7,150) € VSNL) eJ=lel00D 
READ(C7,150) ¢€ VI¥¢s) .U=1,100) a 
95 READIT, 150) 1 Wo¥ts) Jz iio) Sn I ee J 
96 READ¢7,150) @ W3¥¢.) eJ=1,100) oo 
Fs BERD CT TSO) 1 WAV Tl) ET Toor J 
98 READG7.150) ¢€ VSY¥¢J) J=2,100) 
99 READ(7.150) ¢XOFF (uP sJZ=101 NPD 
READC72350) CYOFF CL) JZ102,NP? J 
TO REDD ET, 150) © VIN) JZ1OL NPD 
102 READI7,1SO) € VZ2XCudeSZ101 NP? 
Xt sd= oNP ; 


RRO TT FSO F TWIST 6 EOD QB 
108 READI?7.,.350) ¢ VaXded eJ-IO1 NP) 
105 REDO S00 0 VSR EGY DELO DMB 
READ(7,1503 ( VIV ts) »J=101 NP) a 
io? — READCTSISOD OC V2VCA9 ed=1O1 NP) SSeS 
1086 READ(7,150) ¢ V3V40) .J=101 NP) 
GRADU, S00 vav tab sdz1O1 NPY 


e ee °e@ 
e ~ 
@ a 


ee °® 
Pd 
o 
o 


6 
oe 
© 
oO 





@ - 
210 READ(7.150) ¢€ VSYtudsJZI01eNP) 
TTS WRITEC C1607 TITLE 
@ 112, WRITE (65170) : ee Ne ne ee 
wae are eae IF ¢CICOMP1.E0.0) WRITEC6,180) , 


ss CN ae fee eS oe = as 
e. 


ee ee 




















118 IFUICOMPL.EQ e121) WRITE(65175) 
a 345 IF CIHUB.E0.0) WRITE(6,189) 
~ 386 IF CIHUB ONE Od WRITEC6,190) 
127 (ELN = 
& 116 AL = ALFA = 
9 A N 
120 C THE FOLLOWING SCALES IME OATA BY THE VALUE OF ELNO. 
@ 121 DO 10 I=1,NT ms 
122 RONGIPIXONGI] FZELNE 
123 YONCI)=YON N 
& 128 10 CONTINUE a 
125 00 20 J=],NP 
126 XOFF CID=XOFF CID /ELAO 
& 127 VOFF Cy D=VOFF CI} /ELAD *~ 
128 20 CONTINUE 
ae 129 RETURN 
& 130 C S8eeeFORMATSe¢202 ~ 
13) 100 FORMATE3A6) 
132 101 FORMATOIH »3A69 
® 133 ‘A1G FORMAT(9I49) ~ 
138 111 FORMAT@CIH .918) 2 es 7 ae ee ws tai ee ee per 
_ 135 120 FORMATIOP10F 8.0) ate 
& 136 121 FORMATtIW ,OP1I0FS.3) _ =e “7 ts $< E. so 
ae 137 130 FORMATCOPG6FIO.D) 
138 131 FORMATOIH .OP6F10.2) 
© 139 180 FORMATIOP3F 10.0) ~ 
Tar a ORMAYGIN , Fio.2) Sr ee eS Fits ae tae 
14) 150 FORMAT(OP6E1 3.8) oS 
@  #£Y8f —“Y6U0 FORPATUTHY 80% Se mm —-O———— ay 
193 170 FORMAT(///,2%,°2-C COMBINATION SCLUTION®) 
ag 175 FORMAT(/,6X%,— “INCOMPRESSIBLE WERSION’) 
@ 145 160 FORMATI/ 6X, “COPFRESSIBLE VERS ION®) o 
TRS — £a9 Nee ata La. OF THE FOLLOWING BASTC SOt a - 
147 "Le UNIFORM AXIAL ¢/e9X%4°2e UNIFORM CaoceeLous 72%. 
oe Tas " oe —— —F¥L_ VORTYCIVY ABOUY BOOY I°,7,90, Re WORTICITY ABOUT FO — 
149 3 *pooY 3°,7) 
150 190  PORMNI Ae Nee One Tent ion OF THE FOLLOWING BASIC SOLUTIONS *./,9X, 
@ 15) "he UNIFORM AXIAL % s/,9X —¢°Z2e UNIFORM CROSSFLOW, /,9X, 2 
156e,0—0rtCSsti‘™S 3 73. VORTICITY anouT BODY 1°,7,9X%_.°%. VORTICITY ABOUT * 
153 as "BODY 2%o4e9X%e "Ss WORTICITY ABOUT BODY 3°,/) 
ana ©) elas WEN Ee a te ergy ee SUN phe pec aes ete eet tS ee ge ete gee ge ee eet PY eres ST ee y ot NY” J 
rr a ee he alee ee se a ee Ls ) 
@ re) 
6 a : iS 
& 2 
4 
@ o{_]- 




















~~ 
yy fh he ee So Se ee ey > eS ee ae ~— 
57 VSONCC = ATOTAL/SCRTI2.2) 
& ss IFtIHUB.EO,0) GO TC J6 eee ae Be, eee es = 
59 ASI = yCul-YCL1 
eee 60, AS2 = 2-Y — s 
* 61 IF CC TWeEO 0) cANDe (HDOTCEQ.0.09) WDOTC = &DOTCL*WDCIC2 * 
62 F w 0 .C AND aDOT De Oe 0 wOo = nC - 
63 IFUCIWeEOcOP CAND e FH WDOTCZeFQe060)) BOOTC2 = WOOTC-WECIC) 
eB 64  _—_—s_——s 1G TF UTWeEO-2) GO TO 480 oe 
65 GO TO 50 
66  _—_—_—_s— AO VC = ATOTALSXMC/SCHTU1.0°xXMC 927.0? 
& 6? IF (CIHUB-EQ.0) GO 10 25 eg ee Ver eM ee 
68 VS} = ATOTAL®xMCI/ SORT eXNCIS8 
69 VS2 = ATOTAL@XMC2/SORT(1.04XMNC2902/5 20) 
® 70_ 60 TO 25 ~ 
71 SO IF (IWeEQe2) GO TO 25 
72, VIC = WOOTC/ SGSRHOIOTSACIOCI2. 
eS 73 CALL VBARITCVIC ,ATOTAL -RHOTOT RHCC) ~ 
74 VC = WOOTC/UEPRHOC 4ACI412.0 
75 IF(THUBSEQ.0) GO TO 25 
& ae VICI = WDOTCI/IGORFOTOTSAS}IC12.£ ~ 
17 VICZ = WOOTCZ/IGOREOTOTOASZ)@120C 
78 CALL VBARITEVICI ,A TOTAL ,RHOTOT,R*+0C1? 
e vw ~ CALL VBARITCVICZ ATOTAL »RHOTOT.REOC2D s 
80 VS1 = WOOTC?/1G®RHCC 19451 9912.0 
82 VS2 = WOOTC7/IG*RECC20AS2 9812.0 
& 82 25 LFtvCeGTeVSONCC) WRETEC 6,116) VC = 
oGTeVSONCC) Wt = VSONC oo! 
a8 IFITHUBLEQ.0) GO TC 26 = 
O77 7 RT FECVST CETL VSONCC) BRITETE TITY VS2 oe 
86 IF(VS1.GTe¥VSONCCE WS1 = VSONCC 
OCS WS2CGTCVSONCCD BRETE CG RT OD VSS 
.@ 88 IFIVS2.6T.VSONCC) WS2 = WSONCC “9 
ESET TNS CR Ty de KND.GE.8) GO TO 
90 30 IFUCCKND SE Q.-1).0R. IKNDEQ8)) ELAD = YCU 
QC See TFCUKND.EO. TPLOR. MNOEQ.6)) ELAD = YCU-YCL 4 
92 IFUUKNDeEG.s 03.0R.iKNDE0.5)) ELAD = 2.0 
93 AC = AC/ELND 
94 YCU = YCU/ELND 3 
eae ae a ge ee ee OCR PEND en a TS See 
96 MR1 = XRIZELWD 
a rt a Re ee Te ee pe ee ee ee ee Gg ee ge ee 
98 YRI = YRI/ELND 
99 YR2 = WR2/ELND 
100 XTEST = XTEST/ELNC 2 
Rog ee er ks 1THUB.EO.0) GO TC ee 
102 AS1 = ASI/ELND 
@ TOS WS WSS FEC ND Reyes =e, wd 
104 YCUl = YCUI/ELND 
105 YCU2 = YCU2/ELND 
© 106 YCL1l = YCLI/ELND = 
107 YCL2 = YCL2/ELND 
108 XRH = XRH/ELAO 
) 109 YRH = YRH/ELAD jd 
110 KTESTI = XTIESTIZELAD 
111 MTEST2 = XTEST2/EL AD “ 
eo 42 35 RITOT = RZISOSSTTOTAL | - _ _ hae “LS 
113 “~~ """GRHO = GePT/TTOTAL , 








oS A LT rer CL le _— Tt ewe a ~ ee 


a J __SUBROUTINE CONST SS 

























“ge 5 4 COMMON /COORD/ YON (709) oYONE TOO) .XOFF (200) .YOFFE200)eSt700), 
Be. 1 $1790} ,XxTEST oXTE ST] AVEST 2 oVCL sVEMsYCLI,VCL2,YCUI, Za int: a, Set es - 
a 4 2 VCUZo RRL SXPZyXRH, RI, VRZ oYRHeCUTOFD CUTOF 2 CUTOFH, ee Seta hee ea ee 
'\ 
6 COMMON /CONT/ WC o VST oVS2 YMC oXMC J XMCZ2MDOTCeWDOTC) »bDOTC2< 
& 7 Be on TITLE £3) WINE sALF A A.B eC oD ANC APC AC ARC ASC SALI, a pen ee ee ; = 
8 2 AZT A351 AMT AST yA 1298224832042 AS2, VIC VICL oVIC2 
_ 9 —— —. _COMMON /COUT/ NTANS) NH NP 1W NX sKND s ICOMP om ONXHID eNXHIZ oNXHE 3, Shi eed, et ee ett eee pie a ask 
& 10 i NSTZoNST3NSTT NPERISCISIR AKC SO) oT oP a, Sean eT ET 
1} 2 IHUE 
eas 
Cc 
@ 16 C THIS SUBROUTINE CALCULATES THE CONS TANTS USED IN THE PROGRAM “ 
17 Cc 926560 eee 
i8 P3I=3.341592658 
® 19 mi=0 a 
20 N2=0 
21 IFUIHUBLEG.03 GO Te § 
@ 22 lfi¢iu e£0.0) AND. (W01C.E0.0.05) Miz on 
23 IFQCyY of oOF AND. CWDOT ef 020-4 M2: 
24 IFCCIWCEQLTP CAND. CAMC CEC COOO)) MI 
@ 25 IF C(IN.LEQ, = 















27 IF CCIW CEG e2) oANDe EVS10EO00003) MZ=3 
Oo 8S 01 80=617980.0 7 
29 R21568= 1716.06 ae 
30 6 = 32.178 + — 
e 31 PSTAIC = PSTAT ses 
«SS TRUBS TAT NE 0000) HOTT TSTAT oNE.C.0)) 60 10 10 
33 IF CPT.EQ0003 PT = 2226.23 pe Reet ays 
oe”. Fe AL E020. OTAL = 518.65 3 
35 ATOTAL = S9.COSSSCHTETTIOTALD 


- 


36 CatoT = 1.-0-O.2¢(VINF/ATOTAL }oez 


@ 37 PSTATC = PISCATOTO93.5 ie ae 
NE Pe a i 
3 


8 
9 RHOTOT = PT/ZERZ1I568@TTOTAL) wt Bs cee 
oe  i8  TSTAT = Ftoracecatct ¥ 
4} RHOST = PSTATC/(R2 15 688TSTAT) 
a2 PSTAT = P¥-0.5¢RHO TO TOVING OVINE 
@ 43 ASTAT = &89.0098SQRTICTSTAT) 
~—"s  * BO 
T 










J 
Sc ss i 
4s 10 ASTA = S9.O009SSOR TETSTAT) Bie ee 
CER BROS TS PSLRA TRE EG BB TST RT ad 
47 AMINF = VINF /ASTAT 
CAMINF = f, tO .2 9A PINE O42 


$1 PT=PI1C 


oS TOTAL = TSTAtSCAP INE Se 
5 RHOTOT = PTZIR21S6ESTTOTAL 
@ SSS A = YCU-YeE a an Sc a 
re ee VSONIC = ATOTAL/1.1726 





GO—'-™'"''""'"©°--.. 


- 


whee eR ee 














a Fe Nc een act foto 
a Peas 
is CIOZRT = O.SeRHOST SPT 
& __ 11S VFOAT = VINFSATOTA : 
116 VCORT = VC/ATOTAL 
117 VSIOAT = VSI/ATOTAL 
a 118 VS2OAT = VS2/ATOTAL 7 
Nee eee ence er ae I a 
120 CON2 = 1.-0-0.29VF OAT e082 
@ 121 CONIA = 1.0-0.2*¥S 1OAT##2 em 
122 CON2] = 1LeO-Oe29VSZ0AT#82 
123 RSORTF = CON2¢82,5 
& 124 RSORTC = CON1##2,5¢ -” 
125 RSORTL = CON}LI902,5 
126 RSORTZ = CON21002.5 
& 127? TF UCT WoEQo32 -ORet [bo EQe2)) VIC = VC@#RSORIC ~ 
128 IFCCIVLEQ.1) OR CIb.EOL 29) VICI = VSIORSORTI 
129 LF (CCIW eEQe1) ORI TweEQe2)} VICZ2Z = VS29RSORT 
@ 130 TF UCIW ED cL POReLTheEQe 29) WDOTC = VICSGERHOTOTeAC 41200 ~ 
131 IF CCIW ,EO,)) ORC TboeEQ.Z23) MOOTC) = VICLOGSRHOTOT#ASI/12. 
132 IFUCTWeEGoL POR tT beoEQ02)) WOCICZ = VICZSGSRHOTOTS£S2/12.0 
é 133 JFITHUB.EQ.0) GO TC 36 “ 
138 F = MDOTCICNODOTC2 
135 IFCWOOTCEQ.0.0) VIC = F/IGERHOTCT®AC#12.0) 
@ 136 IFIUDOTC.EQe000) CALL VBARITOVIC ATOTAL ,RHOTOT RHOC} “ 
137 IF(WOOTCEQe000) VC = F/IGRHOCOACI* 32.0 ze 
1338 IFIMDOTC.EQ.0e0F VCOAT = VCCATOT AL 
@ 139 IFCUDOTCEQe000) CON] = 160-DeZeUCOATCS2 so soe 
yao  TFMOOTC.E0.0.0) RSORIC = CONI190z2.5 ts 
14) IF(WOCTC.EQ000) wLOTC = F a 
oe  ~ ~—COT AT—”C— FF = TE -wH OTC? 7 
143 IFtWOOTC1.EQ.-0-0) WICI= FIG SRHO TO0T#AS1/12.0) eo 
1548 IF(WOOTC1L.EQ.0.0) CALL VBARTTIVIC1,A TOTAL sRHOTOT.REOCD 
145 IF(WDOTC1 .EQ.0.0) U1 = FIG SRHOCPAS 19012.0 Saat = 
a4 iwDOOTCi. eOe0) *WSIOAT = VSISATOTAL 
147 IF(WDOTCLeEQ.0e0) CONL] = LeO-Oe Z9¥SLOATOO2 
o-- COG” FF CUO OTED.E0.0.0) RSORTI = CON]1¢02.5 ’ 
149 IF ¢WOOTC1L.EQ.0.0) wOOTCI= F 
150 F = wOOTC-wDOTC? 
151 IFIWDOTCZcEQ.000) VICZ= F/IGPRHO TOT#AS2/12.0) J 
¥OC o£ QO .0.0} 5 AK (VIC2, OFAL -RHOTOT -REFCC 
153 IF{WOOTCZ-EQeDeOd WS2 = F/SIGSRHOC OAS 2 9012.0 
oe 158 SS CWOOTCZ.£0.0.0) BWS20AT = VS2/ATOTAL pe i ey: J 
155 IFIWDOTC2.EQ 2000) CON2ZE = 100-00 Z29¥S 20ATO92 
156 IF (WDOTC2-EQ e000) KSORT2 = CONZ14*2.5 
@ 157 IF CWOOTC2 ck 00060) BNOTC2= F J 
— £58 Sf CICOMP].E0.1) WS) - Vici 
159 IF(ICOMP1-EQe2) WS2 = VIC2 
@  ——S TG CG YF UYCOHPT EO.) ve = vic — 
161 VNFOVC = VINF/VC 
Teo IFUKNOALT.28) WOOTC = WOOTCZELND 
@ 163 TFtKNOLT.6) WOOTC1 = WDOTCASELNE w 
16s.  FFIKNO LTe8) WOOTC2 = WOOTC2ZELNE 
165 PSPTC = 1.0-CLIO2ZRTe¥Cee2 
@ 166 PSPTCL = 1.0-C102R 1#VS1982 sa 
167 PSPTC2 = 1.0-C1D02R 19¥S2902 
168 PSPTIF = 21.0-ClO2R Ie VINFSe2 : 
@ 169 PSPTCI = RSORTFSCOR2 J 
«AGB CPTT = RSORTIOCONII 








ng 


- 


She mee ela aertan te 


am =i 


aman ie 


om sable 


S _ ~ 
saan tieeec ee a 
e  @™©|......eren a ees a a a ee 

17) PSPT2 = RSORT29COAZ1 C—O 
& 172 PSPTCC = RSORTCEeCON 


~*~ 
173 XMINF = VFOAT/(CONzO*0.5) 7 
174 IFUXMC2EQ,000) XMC = VCOAT/SICONI 490.5) 2 
@ 17S TFCXMC1.EQ.0.03 XMCl = VSIOAT/EICCN1190065) ~ 
= @ 
177 QINF = PTe41.0-PSPIIF) 
@ 





178 -QCINF = PIC# (De 7T¢VFOATSRZORSORTIF } oe ~ 
179 QC = P¥#t1.0-PSPTC) 


SCC 5 PICO eTevCoaTeezeRsontc) wee 
181 OC] = PT*¢1.0-PSPTC1) = 
182 QC2 = PTetleO-PSPI¢2) 
183 OC1C = PTC8IO.7eV<s 10ATS82ORSORTID 
188 OC2C = PTCHIO.7EVS Z0ATPO2eRSORT2) i mm 
Ts THETC = TYOTAL/518.69 
186 DEL = P1C/2116.23 a 
187 IF CENHUB.LEQ.O? 60 10 37 “a 


& 
188 WOOTIR = VOOTCISSCRTCTHETCO/DEL eee 
189 WOOT2R = WOOTC2®S CET CTHETCI/OEL 

@ 190 37 WDOTCR 





¢ an 




















= wOOTCeSQRTETHETC) JOEL “ 
191 IF¢ICOMP).E€0.0) 60 TO &} 
192 URITE (6,101) 
S 193 WRITECGol LD VO VKMC,OCePSPIC ad 
198 IFCIHUB.EQ.0} 6O TG 38 
195 WRITECGS121) VSEOXKC1 »OC1PSPIC) 
196 WRITECGSL3ED VS2—NPC2,0C2 PSPIC2 se a ae 
TOT S38 WRIVETE, 181) WINE, MINE OTAF PSP YIF———— be VINE -XMINE, ‘ OS Se —j- i a es 
198 GO TO 51 ee, oh tse 
e785 - 6s —- es aa ee Sac aa era ra al - 
200 WRITE(60110) VC oXMCeOC GCC PSPIC PSPICC,RSORTC 
2 IF¢IHUB.EQ.0) 60 To 39 
S 202 WRITE (6.1203 VS1,x%C1,0C1 ,0C 1C.PSPIC1 PSPTI,RSORT} = = ia 
oe 2030 WRITE 169130) W52 HRC 2,002 ,0C2C, PEPTCZ PEPTS RESORTS aon a wee eae op 
208 39 WRITE(6,140) VINE , XMINF oOINF eOCIAF oPSPTIF sPSPICI oR SORTF >, ae 
oe = ~—S—S WRITETE. praca ee te 
206 WRITE (6,155) ALFA» WNFOVCYVSONIC» VSONCCoWOOTCR »¥DOT IF pWDOT2R —_ 
207 WRITE (6,160) 
& 208 WRITE(6,165) te 
= ee oe 6.170) at, AY, ATC,AS oRHOST, WOOTC, BEOTCI. ubO 
210 WRITE(6,175) VWIC.VIC1,¥IC2 _ 
OQ 21 WRI TVETE ci 809 ee ee we 
212 WRITE(6,155) TTOTALSPT PTC sATOTAL sRHOTOT » THETC DEL 
213 WRITE C6.160) 
@ 214 WRITE (65185) a ere =e a _ ww 
Toes 216 WRITE(G,1700 XA1 VAI xR2,VR2,xTE ST, VL YCULEUND 
216 WRITE (6,190) epee oe ; 
oe 7 STF aT €(6,170) KIESTI,VCLU1,VCUL XTESTO.VCL2 »YCU2 ww 
218 WRITE(60210? CUTOF 1,CUTOFH CUTOF: 
219 WRITE(6,16U) 
eo 220 WRITE (6,200) NT SNP NSIT SNH »KND IW NX oS ICOMP1,IHUB ww 
“221. WRITE(6,160) 
_ 222 ALFA = ALFAe@PIOIBC pe ee _ eee ae 
2230 C888 ee ORMATSe08e8 ~~ ee eg eg we! 
as 224 100 FORMATUIHO,27K,°MACH’ pL IMG DYNAMIC PRESSURE*®.13¥,°FRESSURE RATIO®, a es 
225 1 RX. *DENSI TY RATIO® o/ 414K, *VELOCITY 5 7X _ "NO *oONe INC’. BOX, ‘ 
eo | 226 2 “COMP *o1EXs "INC ® BX—e° COMP ys) ere. O. se, ean tit ae de ates fie eee ____*bud® 
lie agg = ee TOL FORMAT (IHO.2 71, *r ACH 11K, DYNAM IC PRESSURE*S13K78 RESSURE RATIO’. —— se 


— rn er ee 
¢F _ — a —_——_ - = — —— es ee, —__.— ee —— — + -— Te eS Ser 8 eee 


- ‘iatmcte tm of sothhe 


8 




























S am, 
A ETT SE ee 
& an 
228 Cel@X ye "VELOCITY 7X e°NO ag 1 7X e * INC 25K y°INC GAD 
® 229 110 “FORMAT $n, CONTROL * 2IPE 200 302% ZIT PE 200 305%» BPEL Oo 39 3X De IPEL0s 35 - 
230 6X ,1PE1G036/) 
2 lit Ceephatees 
& 232 116 FORMATI2X,*°CONTROL STATION VELOCITY = °,1PEL3e%—_° EXCEEDS VSONIC.* ~ 
233 4,* PRE E = DN RC ATION ° 
238 2 °1S NOT AT THE THROAT RESUBMIT WITH A LOWER WC.°) 
@ 235 LL7 FORMATC2X,°CONTROL STATION VELOCITY P ° EF . a 
236 1 efe* PROCEECING WITH VSI = Te RESUBMIT WITM A LOWER * 
237 ° ® 
& 238 118 FORPATC2x,*CONTROL es VELOCITY = *elPEL3.%e° EXCEEDS WSONIC.* = 
239 1 4? PRO EE TH : Ey mM R * 
2480 2 "¥VS26 
& 241 120 POAWATGsx te) UEROOiaa opaee ake *oEPEID > 3e2XeZCIP E1003 SN SIPEIN. 3 -~ 
242 1 eSXMPeLPELO se 3e6NVIPEIOW 3,7 3 
283 121 FORMATO3X,°LOWER® o/e3X— "PASSAGE oJ PEIN. 392K gS IPEIC 399K» IPEI063 > 
@ 248 1 Z2OxNSIPEIO’ 3/9 ~ 
245 130 FORMATUE3X»*UPPER® o4e3Xy°PASSAGE “pt PELO. 3, 2X p2UIPEIC 3.5K pIPEIO“3 
246 1 e3SXPoIPELOc 2e 6X eo IPEIO ce 3e/) 
& 287 B31 FORMATI3X»°UPPER® s/o 3X%e PASSAGE "oT PEIN 3e2N pIPELC eT 99K ep IPEIOe3s mn 
248 1 2OX,1PE1003./) 
249 ASO FORMATU3X,*FREE % 9/5 3%, *STRE AM “eIPELO 392k o ZUIPEIGe3 SX ,IPEIO.3 
@ 250 1 o 3X), 1PEL10.3,6%,1PE1063,/) “ 
251 191 FORMATU3N» *FREE* ss o3X_ "STREAM * eI PE 1003 92% — IPE 1003 pINe IPEL0 03, 
252 i 20Xe1PE1003,/) 
& 253 ASQ FORMAT(/.9%y PALPHA %99Xe “VINE /¥C° 9 7X y “VSONIC » BX— "VSONICC' IX, Si Ae eee ne, 
258 1 "WDOTCR® 8X oe  WOOTLCR® e8Xe “WOOTUCR®) 
255 15S FORMATI7X,701PE1003,9%)) = z 3 _ 
& 256 160 PORMATALS3 Mga" sienna aeeince sehen e==ea- ease wet ce cee mene saan - : : 
257 Dn te et nen tenn ee ene oe ee ee e- * 7) a 
258 165 FORMATIS,9X%— °TSTAT Xe PSTAT “oOo PSTATC “08K ¢  ASTAT*o9X, *RHOSTAT® 
& 259 1 o7Xe°WOOTC® IX, *MOOTL °o 7X. WOOTU®? ’ 
“360° “ST FO FORMATO 81 lPei0.:.4%))  ””»”»©}©}™©™©™~ ss 
261 LTS FORMATES IN, VIC® g LUX SVICL® gp 10K oe VICU oS oI a SUIPE 2003 ,9KN)) _ 
@ O FORMAT(s.9X, °TTOT® *PTOT* 10% eo °PTOTC® 9X — °ATOT “ohOX, *RHOTOT® » 4 
263 1 BX, °THET® CX, "DEL ) 
268 185 FORMATE TCH oe *XRID LIN VRID oh CX * MRIZ® IN, PYRIS  oSKy  XTEST * ol OX 
® 265 1 eo YCL*shIN, *YCU® LOK, *LND®) é 
56 7U DM oIX, °X e1OX, °YCLI 0 LOX, *¥CU e8X,°x = e1lOX, CL 9 
267 1 10X,"ycu2*?) 
& 200 FORMATU/, 10% 4 °NT® 5 EX» NP o5X_°NS1°,5N DWH o5Ne MND, Ske luc 5K, ud 
269 i TNX SK o*ICOMP1® 5X o* THUB % 9/5 SX pT3 yO ol 3p SX pl 3p 8X op 139 5Xo 
270 2 TSH oLTeSNelZoTNolboGXel ID 


271 210 FORMAT¢/ > 7Xo°P-S CLIOFF 2° 3X e°P-S CUTOFF | 2 P-S CUTOFF 3°,/, 3 
BX ,1PE10.03.5% -1PE10.355X%-eIPEIO.3? 
273 END 





| 
| 
| 


® 
U- 








je 


wtade. 


-e~ 


we. 





No SU Oe et 
2 COMMON /COORD/ XON 4700) ,YONC 700) .XOFF (200 D,YOFF 4200) .S¢ 700) 
3 ______1__ SIU 700) sXTEST XTESTIgKTESTZ29VCLeVCUSVCLIeVELZ2sVCULe 
& a 2 VCUZ XR oXRZ2eXRH > YR1 eg VR2 pVRHCUTOF 1 »CUTOF 2 CUTOFH, = 

















6 COMMON /COUT/ NT,NS1 oNH oNP oT WeNX SKND -ICOMP oK ONXHI1 NXHI2ZoNXNE 3, 
& oe 1 NSTZoNSTSoNST7oNPFRESOD TRAKG3OD aM) MZ eICOMP] oIPL i 
B 2 IWUB 
9 COMMON /CONT/ VC oVSE pVSZeXMC gUMC 2 XMC2eUDOTCANDOTC Le bOOTC2, 
@ 10 1 TETLE C3) gVINF CALF AA eB CoD pAIC cAZC eAIC ARC LASC CALI, ~ 
AZLoA ASLA a A32,A92,a85 
12 COMMON /CONDITS TTICTAL PT sPSTATSTSTAT PSTATC, ATOTAL oPTCeRHOSTs 
& _ as i RHO TOT AASTATLQCINE 
Fg OSE aaa sana cw wenn nnc wan ween eccnuceae suaee CNT DOO en Oe Ee OO 
Siesta Do c 
& 16 C THIS SUBROUTINE FINDS THE HIGHLIGHT ON EACH BOOY AND ~ 
17 ALCULA AREA OR A ON- A OF F-BODY POIRK 
18 c 
eo 19 ; PI = 36191592658 “ 
20 Cc 
— 21  _  € —=—ON-BOOY POINT CALCULATIONS 
BE a wh pe ee eae ee ee eo “~ 
23 JJ = NS141 
28 JJJ = WH * 1 
e 25 NXHII=1 “ 
~ 26 ~ 00 16 J = 2eNS) 
27 IF UXON GJ) GE eXON(J-13) GO TO 10 3 
eq —_ % = ge IN a ee ee a 
29 10 CONTINUE 
~——s« SO”*~“‘=‘; S.. ©6cSS EOD CO TC 2) 
& 31 J1 = JJ + 1 ey 
32 ue = Jde-l* SO OS SSS SS SS SS Se —— 
33 00 20 J =JleJd2 Sea 
Qe ae AF AROC 0 sGe MONG“) 60 WO 20. SS oo 
35 NXHI2 = J 
oe 36 a KHIZ=XONGNXHI2) 
$ 37 20 CONTINUE ug 
38 pois 2t 32 = JuJJei 
39 IFC IHUB.EQ.0) XH12-299999.0 
oC ee ee a 
41 IF CXONGJS) oGE XONIU-13) GO TO 30 
oF &2 NXHI3 = J 
®@ a3 30 CONTINUE Ww 
NS TSENKHT YO as " 
a5 CALL SURF 
os ws 
a7 C OFF-BOOY POINT CALCULATIONS 
a8 Cc 
@ a9 K =0 Y 
—. $0 NPI-NP-1 
51 ; 00 105 I = 1,NP1 
& 52.C~<C~‘Ct;«‘t;‘  FRCXOFF(T).UT.xHI2) 60 16 101 ~ 
53 IFGI.EQe1) GO TO 101 a 
58 L=Le) é 
@ 55 IF(LsEGoL1) GO TO 102 _ ; aoe ee Se ee 
56 AA=VOFFII*1L)-YOrFF tI) = : 


See = a ee a er ee . a aes eee ——— 
ns a 
ae eS re re ee 


ale 


te, 


iE a 


~ oe Se 


Petre atet 


oe 


—aletieliie ej Rmseb 














& -~ 

& a a aa a a 
57 ABZYOFFUI)-YCFFtI-1? 

# 58 IF(ARSCAA-ABD.GT.0.001) GO TO 102 : Se Nays, ~ 

OFF I1) eFC NOFFile1)) GO 10 eS eee 

59 101 IFUXOFF UL) .EC.XOFFIIe1)) GO TO 1C® 
60 102 K=Ke) 

& 61 NPPRIMD = 2 oe 
62 L=NPPREK) 
63 tiztel 


& 6% __105 CONTINUE Seeks Mes ne 
65 KoKe) eee Rt Be ne SS nr ee 
66 CO 
& 67 DO 130 I=1,NP - 
68 Ic =0 
69 IF¢CIoLENPPROLD) M421 
70 PEATOLESNPPRELDD GC TO a2 le 
71 DO 111 J=2eK 


® 1$ TFUIHUBONE.O) CALL INTERERON »VOM gNNONXHI2 XOFF II), 153 “~ 
? IFCXOFF¢1).6TeXHI23 60 TO 180 
78 IRAMUKRIZ 3 
= ~~ F9 60 TO 110 “~ 
O IFtYOFFtI dt T.YS) €0 TO 190 
IRAKEGK1) = 8 
is 60 TO 110 7 7 
90 FRANGKID = § 3 
60 TO 110 pac 
€00 IRAMCK]1) = } 
8 D CONTINUE 
ay RETURN 


@ 88 END £5 Sd cage 
I 




















ae TP ee ee hae re ene peer a a a eee ras ee a 
ns A a a 
® not SUBR NE AN i ee eh al See ee ea ee -; 
; CCHMON /COORD/ XON (7U0) pYON( 700) sYOFF (200) oVOFFO2OCbeStTO0)» 
os A SLEWOD sXTIESTeXTESTI XTEST2Z4¥ GhavGmaVCLLaVCb2eVCUle 

& YOU ZX 1g XPZyXPHe VE] o YRZ oVRMoCUTOF | g€L TOF 2 oe CUTOFH = 
6 COMMON /COUT/ NToNSI NH WNP, IW NX END oICOMP GK WME S CR DMI2 NANI 3 5 

& I } ————NSTisNSTSeMSTTSNPPPUSCIeIRARCSODaMRaPZeICOMPLEIPL, mets aes 
8 2 THUE 
9 ¢ fn eee @ a @ Cease nweee eee es 2 = oe wae Sete O28 Bee es 22S SSD TEBE BOWes 6 & eS @2@oe ea @ @P@ ee ese = = 

& 10 C Cg 

cy yw a ri ARS id 4 > ed : 

12 c 

eo 13 NS = NSi-] - 
18 DO 10 I = LNs 

15 IFUXONET*O2RDCEQGANONITI? ANGET I= 

@ 16 IFEKONCIO12cEQeXONGI)) GO TO 10 ” 
17 ANG{I) = ATANC(YON IZ 01) -VONET)D/ EXON C202) -KON ETD) ) 
. 10 CONTINUE 

e 


NI2 = WH-} a eT 
20 IF (1HUB-EO.0) 60 TC 21 
os ae N = MSO 
33 00 20 I NeNT2 “~ 
23 IF(XONII01) EQ.XON ITD? ANGUT3=3-19159/22.0 7 
28 IFOCXONUIO1).EO.XONIID) GO TO 20 


@ 25 ANG(1) = ATANUCYON (162) -VYONET DDS UXON CIOL D-XONEIDD) a re 
TT CT wy 
27 22 N = NT202 os . 
® I A ie 
29 DO 30 I =N,NT3 weet. he 
+ e e ( =32]4 Joe 
r 31 IFUXONUT+1)0EQ.X0ON(T)) GO TO 30 1 


SE ERC TTY = AVANT TVOR TY OT) =VONT TVD AXON CT OTD RON TTY SSeS 
33 30 CONTINUE eth ee, 
eo 3s URN a 
35 END 
A 
e 





eo | oe eee nc 
@ oO Se ge ee ee ee ee ee ee 
@ ees Pl 
& Bg ge ee eS ee i a ae a 2 
@ oo ey old) _ 


























poe Deca ROL Nand a Soo, eet: Soh tee Be cet Nias, oe wis Mao Nis nd Pe ene ied 
one ae aes SUBROUTINE SURF _ es 7 
2 CORMON s/CUORDsS XONI7TO0) -YONE 700) .XOFF ¢200),YOFFL 2003 eSt700), 
3 1 $21¢200)  XTEST YTE STI pXTEST2 ,YCL »YCU ,YCL J, YCL2,YCUI 
[ a 2 YCUZe XR] MARZ eXRHe YF) 5 YRZ ep YRHeCUTOF! ,CUTOF2,CUTOFH, = 
5 3 ACA ANG R R 
6 COMMON /COUT/ NT,NSL NH NP IW NX KND pS ICOMP KK ONXHI? NXHI2,NXNI3, 
> 7 1 7 NSTzoNST3yNST7 pNPFEL3C Dd pIRAKUSOD SMI M2 gICOMP] IPL » = Z 
8 2 IHLE 
_s ee COMMON/FINDERS JCHK) pJCHM 2 pJICHK 3 
ey 10 Crt nn nn nn 2 ne = oe oe oe ee - wow mew oe nw 5 en oe ee ee = 
12 Cc 
12 C THIS SUBROUTINE CALCULATES THE SURFACE OISTANCE ON EACHK BODY 
oO 13 C FROM 11S HIGHLIGHT sien 
18 Cc 
15 CALL CHECK() NS ,TRT,XOND 
oO 16 ZRS=XARI = 
i? IF CXONGIRT Pe GTeoXR1} NRIZXONEIRID 
18 DO 2 I=2,NS} 
19 TF CCC VONET) GEC VRID ANDe CYONEI— 1D LE pVR}PDoORALIVCHATD OLE VARA De = 
20 1 ANDeCVONCI“1LI-LGEIRIDI De AND CIC XONE I 2eLE oXR1 DAND e CXONCI~-13.6E- 
21 2 XR1IPP OR SAC XONET ) GE oXRI DANCE NONE I-“DDCLECXRADDDD JCHMI = I 
Ee 22 2 CONTINUE - 
23 BRI=ARS 
28 IFIIHUB.EO.0) GO TC 5 
oi 25 KKONS1O2 = 
26 R2ioKK—-1 w 
2? CALL CHECKIN 2,NH,I&,XON? ___ os — 
CO 28 XRS=KRH : 
29 IFAUXONGIR DG TeXRH) XRHIXONIIR?) 
30 DO & I=KK,NH 
- 31 IFCUCUL YON (I) GE eo YRED CAND e (YON ET — JD eLE YRHDPOORGUIV CHAI D OLE oYRHD 0 
AND. CYONCI—-1 b06E€ « YRH)) Do AND es 000 BONE I DeLE oXRH) AND « CXONCI-1) -GEo 
33 2 KRHDP ORC EXONEI GE eXRHDANDoUXONE I-12 Pe LEXRHDDD} JCHK2 = 1 
<= ONTINU 
35 XRHIKRS _ 
36 5 NKINHe1 
C) 37 MK 2=NKO] - 
We etre oN RR a CALL CHECKINK NT, IR2,XOND 
39 MRSIERZ 
Cc - & (TR2.F0.0) = . 
41 IF(XONGIR20 eGTeNRZ) ARZ=XONCIRZ) 
a2 00 6 I=NK2,_NT 
a3 IF CECE YON UL) GE oe VR Z) cANDe (VON ET ~29.LE oVR2ZIPOOR ELV CACTI OLE oYR2 Do 7 
TO TTR o TYON IT o UCUXONCID.LE oXR2 POAND «UXONCI-1)0GEn 
a5 2 ARZPP COR SHE XONEL) GE eXRZFCANDe CXOCNET “12 0LEoXR29993 JSCHKS = I oie sais 
aT ARZ=IURS 
as IF CYONCUCHKI ).GE-YR1) GO TO lO 
: a9 SCUCHKAD= -SORTIIXCNEJCHM) D-xX RL} 492% (YON EC JCHK 1 D-YR 19982) _ oo : 
Tae eS Vee ot F NNi - JCHKi¢2 fect San = Poa 
51 = NN2 = JCHKI-1 noe ete nes See Atoka et ee 
= as ~SENND-2P= SORTCUXCN(NNI-1)-XR1 D820 CYONENN] -2 )-YR ] 3992) 
53 GO TO 20 
54 10 SEJCHKID = SORTCEXCNEUCHM1TI=-XR1) 0926 LYON ET UCHK IT P-YR 11082) 8 
oo NN1 = JCHK1¢1 hosel te hee ae Son. eae seas en Pe Te i 
56 NN2 = JCHK1-2 i 
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Mi E20 OO SO N= NNDGNSD ere 
5¢ StI) = Sél-2 D+SORTCIXONCT D-KXONCT-1)3982¢CYONE 2 b~VOR( 1-1) 3) 902) 
OD 30 CONTINUE geste 4c ee ss 
& 61 DO 4G II = 1,NAN2 
pst ce+l- 
63 StL PD = SELODMHSORTECXONET D-XONCT 619) 002 9¢ YONA I I-VORCTO1) } 002) 
ie OW ct te NG CONTINUE a oo oe ty 
65 IF cIHUe.EC.M} GO TC 81 
66 CF AYONECICHK2D.GToYRH) GO TO 50 
& 67 StJCHK2) =-SORTUUXCNUICHK2)-RH} 9824 LYON ICH 2 D-YRE POe2) 
68 NAL = JCHK 202 
69 NN2 = JCHK2-) 
x 70 SUNNI-1) = SORTLAXCNENNI<13-XRH) 0920 1 YON UNNI~1 D-YRE 99929 
71 60 TO 60 
72 50 S4JCHK2) = SORTECXCN CICHK 23-XRH 10924 LYON ( JCHK 2)-YRE D082) 
@ 73 Nh1 = JCHN2Z¢61 
7a NN2 = JCHK2- 
75 SINNZ¢1) =-SORTOCXCN INN 2¢23-XRH1D092¢ CYONUNNZ41D-YRE 1002) 
® ae NOE. 60 IS = NSIe} 
77 1$2 = NH 
8) DO 70 I = NNW1,IS2 
& 79 SCL) = S¢I-1)¢SORTAIKONCL D-XONET-1) POR2>CVONECTIP-YORCI-1) 2002) 
80 70 CONTINUE 
Bi DO 80 It = YS yNN2 
& 82 I = WN2°1$-I1 
TEE eS SUV = STI+IF-SURVITKON P-KRONCL Ol) SSS OT VOMUY S-VORUTOL VTEC eee 
88 80 CONTINDE — 
eq 55 BEET VON TITRA T.Va27 60 TO 90 
86 SCUCHK3) =-SORTEC CNC UCHK 3 D-XRZ 30020 LYON CUCHK 3D-YR 23902) 
8 WHI -— JCHK3¢] nN 
88 NN2 = JCHR3-2 ° 2 
Te ad ee  STRNZeTF = SORTTIVCN INN +] I-XR2 30020 (VON CANS ¢ 1 D-YR 2007) = 
90 60 TO 100 mo es 
oe ———Ty_ $8 StICHKSY = SORTUUXON CICHM US—KR OV SOS TVONM ESCH SU VR UOC ID ——< 
92 WWI = JCHK3¢2 eas 
93 WN2 = JCHKS3-1 
& gee SENN1-2) =-SORTCCXCNENNI-2 D~XR2) 0826 (VON CWNI-1 -YR 21982) wae 
2 —_ 7 oy EO ee RAL ee ee ee te ig er ce wee de web See ee ema ey 
96 DO 110 I=NN1 NT 
ns > = =TI-50 
98 110 CONTINUE 
99 00 120 I2 = IS,NN2 
& 100 I = WN20IS-T1 
= OE “> eta = . 
102 120 CONTINUE 
eR a ee ee eee ee 
104 END 
Ors 2 see eee 
@ 
7 As =: _____ SUBROUTINE CHECKUIS, IF ,1IRT xO) = a 
2 DIMENSION xX01500) 
ae cee et 1¥=0 
_ ® 1ST=1501 
5 00 5 I-IST,IF 
6 IFCIT2EQe1) GO TO §& 
OTR EMOCTDGTeKOCT <1) IRT=I“2 ee 
8 IFCXOCIPCToXCUT-1)) It=1 
seby cette 2 oP ie ye US CE ONTINGEL a ou 7 es Sok 
@ 10 RETURN 
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1 S1t700) »xTEST aXTESTIAXTEST2 eYCL eo VCU aVCL Le VCL2eYCUR, 
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COMMON /VELOC/ ¥11700),.¥21700),.V2170U).V4t7TN0) VS4 ICC) V1IXI200), 
i st VZ¥I200) pV3XE 200) ,vaxt200) 2V5X(200) .VIVI200}, ~ 
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COMMON /CONT/S VC VS! ,VS2_XxMC pXMC 2 
10 1 TITLE C3) oWINF oc ALF AA eB CoD AIC gA2Ce F3CSASCSASCAL Lo 
AS 
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23 PSOPT (700) .CMHACHE 300) »XMACHE 100).CP 11700) ecPCi700), 
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27 C THIS SUBROUTINE SOLVES FOR THE COEFFICIENTS &,.8,C, ANC C AND Ses - 5. ae ae a 
@ 28 C SOLVES FOR VBAR ANO V INCOMPRESSIBLE FOR BOTH ON- AND OFF-80bY ‘ >= : = a 


29 C POINTS ps 
See Gy ee a —_—~ 
@ 31 OIMENSION YWEWI50) VNEW 1650? .VNEUZISODSYNEWSESO) .VAEWSESOD, a 
See ee i ~ WNEUSESC) ; " a . pee ote oe ee 


33 6 = 32.178 
CN FF E001) G0 TO 30 ta 


35 C__THE FOLLOWING CALCULATIONS ARE FOR THE CONTROL STATION 
— «S6”””~”~CS”S*C* SPREE OF THE HUB ee ee 
@ 37 A = XTEST rw 
TI = 


oe 3 7=CS«ws 0 16=Nx N51 — @® 
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6 CONTINUE 
@ —— ALL IN 2UNMHIT OAS 1 oA OVE VIL DV CL OV ILSVSLOVSLD w 


87 CALL INTERZ CNHI NXT 3A oVH VIHVcHeV IH, VaH,VSH) 
“8 00 10 IF = 2K 


@ _ 49 IFA(IIEQe1) I=] ~ 
50 IFdIITCEQ22) GO TO 33 





7 51 I=NPPRETI-1) 4) a 
52 ER TETAS NE CKOFF AID) CO 10 10” “a t ww 
, 53 IF¢TIILEQei) NPL = } 
54 IFCIT.«EQ.1) KPH = APPREDD 6 
@ 5S IFCIT.EQs1) U3 = NFPROEDD he 
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a: ge Atte ee eats = Bee er a oe ete ee 
5? IFCITeEGes) NPL = 1 ap Ee, Pe ee a ES ee 
= ~ s8 TF tLTsEQ.J) NPH = NPPREJI-NPPROS-DD 
59 IFCIIT.EQ.U) JS = APPREYD 
60 SS CONTINUE ee ea te : 
a 61 10 CONTINUE 
62 IF (NPL EQ.) WRITE (6,1000) : 
63 IF (NPL eFO.0) STOP 2 
@ oe 0 0 I NPS BN hw ee Ge ee ale be 
65 YNEW(JU*2-NPL) = YOFF Odd? 
— BC EWLtUJ*2-NPLD = WAXES) peek at Rte to ee 
s 6? VNEW20Us*2-NPLD = W2xluu? ey 


eyee ot 


68 WVNEWSIUJe2-NP x 
69 VMEWSE UI 2-NPL D WAX tld) 


® 70 LFITHUB NEO) VNEWS(JJ*2-NPLD = WORE yD 
7 20 CONTINUE 
72 YNEWED) = YL eee 
® 73 VNEW}CI) = VILO@COSI(THEALD oO 
NCS TD 
75 VNEW3CL) = V3LOCOSITHEAL) 


ee i VNEWSEL) = VaLeCoS (THEAL) a 
77 pranene 


IF¢(THUB NEO) VNEWSCL? = VSL SCOS (THE ALD 











_ 78 YNEWINPHO2) = YH eee het 
OT WE ML ONPH O27 = —WLPOCOSCTHEAMD “ 
———-____§0.. WV EZ IPH) = @V2POCOSETHEAHD — 

al VNEWSENPHOQ2) = -VIFSCOSTTHEAH) 
e 82 Baca UN CE ie De WE COS UE HU tao a ee eh RO Me ee as 
— es IF¢CIHUBeNE Dd VNEWSENPHO2) = ~¥5F8COSITHE AH) we 
88 WPT = NPHe2 sie. = re 
eG 85 ERLE ENTER CVNEWT. YAEUAAIC NPT) . 
96 CALL INTEG(VNEWZ, VRE ,AAZC NPT) 
& 88 CALL INTEG(VNE WS» TREN SAASC NPT) 2 ae 
89 TF CIHUB NEO) CALL INTEGUVNE WS VRE We AASCONPUD 00 
90 AROFFUNPL )=YCU-YCL Seeeeieeas! 
Qo si AICS ASICCRROFEINSL YY °° CC a ee ee } 
92 A2C = AAZC/AROFFERPL) Ss te 
93 A3C = AASC/AROFFOCNFLD 
@ 98 ASC = AASC/AROFFUNFLD : Pennie 
BTR CTHUG NE OD ASC = ARSCTRROFRONEUD 
96 IF(IHUB.E0.0) GO TC 215 se a ae es 
oe 7 7—— «TE FRE FOLLOWING CALCULATIONS ARE FoR inf Low CONTROL <TATION ae ‘3 
98 30 IF(M1.E0.1) 60 TO 40 
99 IF(MZeEQ61) 60 TO 70 
@ 100 80 A = XTESTI = _ ll 
i. 101 ~~ NPE = O 2 
102 NH1 = NSLo] an Dette ee ee ee eee 
@ 105° CALE INTERZONST SON) oh YNEWCLDSWAFU Ci) WNEW2 C1). VAER3UL) eo WNENSED) teat 
108 1 oVNEWS (19? ae 
105 OO 41 TLO=NST3.NS1 
@ 106 IFCUASGEeXONI TIO} ANCOUACLToONONITIOCID DD THEAL = #ACCTIOD Mee ie 
107 ~ @T CONTINUE 
a ROB he MOONS ROSNU LINN G 2 one ie ey be gah oe re a 
@ 109 TFCCASLEs XONCTIO)) ANOS’ (AcGTeXON (1100193) THEAH = ANGIIIOD ~ 
110_ ‘42 CONTINUE 
PAD WNEWTCD DY = VNEW DODO COSUTHERE TTS 
@ lie VNEWZEL) = VREWZ2CL@COSETHEAL) Se ees a ee Ses pe es, ot eee et 
~NEW TOT ISCOSTTHEALD 
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& 115 ies VNEWS(1) = VNEWSCLI®COSCTHEAL? Ba a a Se 
116 CALL INTERZCNHI gNXEI 2 0A ce YH eV 1HeVchy¥ SHeVEEVSH) 
117 OO 50 IT = 34h oe pei Tee kt 
e 118 IF(IY.EQ0)) I=] 7 
1}9 IF¢LLeEG.1) 60 TO 45 
120 I = NPPRITI<1)01 
@ i123. __95_IFCIAsNEXOFFIIID CRC YOFFEI).GT.YCUI1)) 60 TO 50 7 c >, foes es sae 
122 IFCITEG21) NPL = 1 
123, __sC FAT EQ,1) NPH = APPRODD Beet 
@ 124 IFCII.EQ.1)) «JS = APPREDD ™ 
125 DO 86 y=2 0K Saas neeeemnemnrneeeee eee 
126 IFUTEeEQ.J3 NPL = ARPPREJ-1)01 
@ ev. 5s. IF(TEsEQeJ) NPH = NPPREJI-NPPRIJ-1) -~ 
128 IFCITe£0.0) J3 = APPREDD 
129 86 CONTINUE sD neereeemenemeee ee 
@ 130 SO CONTINUE = 
131 IFCNPL EQeO} WRITE 16,1002? 
132 IFINPL.EQeOF STOP 
® 133 = DO 60 JJ = NPL.J3 a 
—— 2. YNEWCJIO2-NPLD = YCFFEUUD 
135 VNEWLIUIO2-NPLD = WEN CUS) 
ee 1360 WNEW2 CSO S-NPLD = W2xtduD = Se a ee ae oe 
137 VNEWSEIJO2-“NPLD = WINE UU) a 
138 VNEWS(JJO2-NPL) = WaXtdu) 
€ 139 VNEWSCUSO2-NPL DD = WSXCUe) ue 
7 = Are oo ot SR ee ee is 
181 YNEWENPHO2) = YH i 
OQ NR WREDT TWP SZ = VIFSCOSTIRERRY Oo - — a ooo, ee 
143 VNEWZUINPHO2) = -V2ESCOS(THEAH) 
~ 14g WNEWSENPH*2) = -VIFSCOSITHEAH) 
2S 195 VNEWSINPHO2) = -VarOCOSITHEAH) eet 
SO TRG VHEWSTRER ZT = -VSPOCOSTIHEAHD 
147 NPTINPHO2 
rr Tae. ae — — CCALU YNVEGUVNEUY. VREW, ATI NPT) te 
149 CALL INTEGIVNEW2,YAEW,AAZIL NPT? 
150 CALL INTEG(VNEWS »VAEW,AASE NPT) 
6 151 CALL INTEGUVNE WS, YAEW,AAST NPT) s 
TS CRT TIN VEG CWE WS 9 VRE Wg ANS Tg NB = 
153 AROFFOUNPL D=YCU1-YCL} 
O77 TE STA ROFE TRL ) 
155 A211 = AAZI/M@ROFFUNFL) 
156 AS1 = AASI/AFOFFUNFLD 
eS 157 AY) = AASI/AROFFUNFL) wd 
~~ USB AST = AAST/AROFFUNFLI ” 
159 C THE FOLLOWING CALCULATIONS ARE FOR IFE UPPER CONTROL STATION 
@ = ~— «18 TF UTHT E020) AND. th? 06020?) GO TC ab wd 
161 A = XTEST2 
162 NPL = O 
& _163 NH = ONHO? tos 
68 oe NSH= NH 
165 CALL INTERZENXHI2 oNSHoA SYNEWEID oS VNEWIELDOVNEWZ ITD WNEWSOID, . 
oe 7 Oe —> : VNEWR EID VNEWS OL DD w 
167 OO 73 110 = NXHI2,ASH 
168 IFCURLGE XONCTIODDAND. MALT OXONI710019)) THEAL = #NECTIC? + 
& 169 71 CONTINUE _ _ = Sai, ca __ pee one eee Ae! e! eo * he * 
TS D0 72 T10=NH1,NNHT:? 
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a a ee 
& = SE ae a a eR Eg Ee ee, a er oat 
1/1 IF GUA.LE .XONCT103) AND. (AGT. XON(T10 2199) THEAH = SMETEIOD 
& 172 72 CONTINUE ~~ 
. 173 VNEULC1) = VNEWIE1 I@COSCTHEAL D 
174 VNEW2C1) = VNEW2E 1 COS (THEAL } 
7 @ 175 VNEWSCLD = VNEW3ET D@COSETHEALD A 
176 VNEUS¢L) = VYWhEWS G2 DPCOSETHEAL) 
177 VNEWS¢1) = VNEWS C1 DOCOSECTHEALD 
@ 178 DO 90 JJ = 36k i 
179 IF¢II.EQ0.1) 121 
180 IFCIILEQ.}? 60 YO $5 
@ 18) I = NPPRELI-1)01 a 
182 95 IFUtA MEXOFF CLD) CRetTRAKITI INES?) 60 10 90 
183 IFCII.EQ22) NPL=1 
eo 188 IF(ETEQ.2) NPH = APPRIYD * 
185 IFCII.EQ.1)9 9 J3 = APPRI1) 
186 DO 96 J=2 xk 
i @ 187 IFCTILEQe J) NPL=I a 
188 IF (IIT EQeJ) NPH = BPPREIJP-NPPREJ-2) 
189 IFCII.EQe J) J3 = &PPREJD 
i; © 190 96 CONTINUE “ 
191 90 CONTINUE 
j 192 IFENPL.EQe0? URITE 16,1002? 
eS 193 IFUNPL«EQ.0) STOP “a 
198 WPT = NPH*2 Zs ts 
; 195 CALL INTERZIWHINXET 3 oA YNEWINPT DOYNEWLENPT?  WNEWZ CONF T),» VNEWSOENPT) 
; 196 1 oVNEWS INP TT? .VNEWS INP TT)? ie a 
ys we VAEWICNPT) = -¥NEW DENPT)@COS ( THE AH) ne 
198 VNEWZINPT) = -YNEWZENPT DOCOS I THE AH) eet ints 
@e “199 VNEWSINPT) = -¥NE WZENPT 20COS ( THE An) “ 
200 VNEWSINPT) = -¥NEW&INPT 29COS (THE AH) a 
201 VREUSUNPT) = -¥NEWSENPT DOCOS ¢ THE AH) 
Ss 202 DO 100 JJ = WPL,U? __™ 
~~ BGS NER TIO S-NPLD = OY J 
208 VNEWL(CSUJO2-NPLD = WIX UU) alee. 
Se 205 NEW? TiJded- = W2xTID) a 
206 VNEW3SCUS*2-NPL) = W3xXtUa) 
. 207 VNEWSEJI*2-NPLD = WOxXtUS) 
1 @ 208 VNEWS(JJ+2-NPLI = WSXfuu) tab 
ee Oo “100 CONTINUE 
210 CALL INTEGCVNEWL, YRNEWSAAL2 NPT) 
oe 21f°° = =—— CALE INTEGOWNEU2 oYREULAAC2 NPT) w 
212 CALL INTEGIVNEUS »YREW pAA3S2 NPT) 
213 CALL INTEGIUVNE WM» VRE W AAS 2 NPT) 
@ 214 CALL INTEGCYNENS »YAEW,AAS2 NPT) or) 
OTS ROP ONPC ID =VCU2-¥CL2 
216 AlZ = AAL2/MROFF ENFL I 
1 QT ST AROFF UNEL ; ww 
218 A32 = AA32/AROFFENFL) 
219 AN2 = AASZ/AROFFEUNFLD 
. & 220 AS2 = AAS2/AROFFUNFL? we 
; 221 @O IF (M1.E0.1) 60 TO 200 
; 222, IFIM2¢E021)? GO To 210 ae 
. @& 223 C THE FOLLOWING IS FOR THE CASE WHERE UPSTREAM AND LOWER we! 
1 224 C CONTROL STATIONS ARE INPUT 
225 VCkL = AYC-ASC A 
. ae WELZ A31-A0) se a te ete Sn ee ee 
; = 227 VUl = aSi-agl , 
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228 VCU = ASC-ASC CE 
@ 229 &_ = VINFSCOSIALFA) ms 
230 = VINFO@SIN (ALFA) 
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231 D= ¢€¥SI*V¥CL-VCe *( ACAICOBOA SVL I -(CACATIOBCAZ LI OY / 
@ 232 1 vCuevL)) mm 
= - * eB -De 
238 60 TO 220 





& 235 C THE FOLLOWING ITS FOR WE CASE WHERE UPPER ANO LOWER ~ 
236 C CONTROL STATIONS ARE INPUT 
237 200 A = VINFSCOStALFa? 
@ 236 B = VINFSOSINGALFA) ~ 
239 Vi2 = A32-a82 
240 VL1 = a3l~aal 
@ 241 Vul = aS1-a91 od 
242 VU2 = aAS2-A82 
243 O = (VSI*V¥L2-¥S2ey 2+tAal2>Bea2200VL1-tA0A11 0B OA2 1 DOVL2 JivulevL2 


248 i ~VU2S¥LI > 
245 = (VSI-GAeAL oB OA -Oey v 


@ 
296 60 TO 220 

@ 2a7 C THE FOLLOWING TS FOR WE CASE WHERE UPSTREAM AND UPPER me 
248 C CONTROL STATIONS ARE INPUT 
289 210 A = VINFe@COStALFa? 





@ 250 ® = VINFSSINGALFA) a 
251 VCL = A3C-aac 
252 VL2 = A32-AG2 

@ 253 VU2 = AS2-a82 ae “a 


255 O = _UVS2*VCL -VCOVL 2¢ (AGA 1C OB CAZC DVL 2-109 A12¢B9R22 DOVCLDSEVUZOVCL- 
® es 1 WCUSVTEZ) = aia 


257 C = 6¥S2-(a8A129802223-DevU2 D/VL2 
— a eee 60 To 220 
@ 259 225 A = VINFOCOS CALF AD 
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260 B = VINFeCSINGALFAD 
261 C = (VC-APAIC-BOA2ZCE/IAIC-A8C) 
@ H 0 N RE THE CALCULATYONS FOR W AND WAR We ay 
263 220 DO 230 I=1,NT 
IF CIHUB.EQ.0) VINC CL DTACVICT DOBSON 2UT DeCedVStT)D-VelTDD 
@ 265 IF«(THUB.EC.0) 6O TO 230 wh 


66 VINCCID = ACVICTD BO V2C TP oCeCV3t1)-Vatl)) eOetUStl)-Val ) 


267 230 CONTINUE 
OB BRT — ® 
269 TFUTHUBCEOO) VXINCUTIZACVINGT) CBOV2NETD OCOLV3XEI)-Vaxel?) 
270 IFCTHUB SECO) VYINCULD=AOVIVG1)*bO¥2 YELP eCOCV3YUL)-VaYEL)) 
e 271 IF ¢CIHUB.EQ.0) 60 TC 235 


__w 
272 VMEINC CLOSE ASVIXCT DB OVZ0TTDeCetv intl )-VaN TID) soe tVERTIDcVexTIII a 
273 
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oe §} 2% —  @S5 WREITY = SOPTIVEINCET Dee2*VVINCClisea) ae Fe Ss el are Oa ee a 
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71 WRITE¢6,521) —_— 


ns ie She ae DO &S & = NPL oJ3 7 
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93 100 CONTINUE 
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"95" ~~ ~~ ~@p wRITE(6.500) I 
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oe ©) SEE ERG UIT PO pe sage te AM, Nees ge ay Ee Se aw ey We ee Oe a ~ 
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oe 325 NN-NH*¢] ~~ =. ee - _ zg 
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4 C AND RETURNS IN CD THE CORPRESFONCINE VALUE FROM Bo 


NIEFFOLATES ARRAY A FOR THE VALUE ¢ 


oi eh a ts et me 


a ee 
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248 W3 = CCHKEMA DPIC CHK AMDIZ EMER SLD KEM H- LIDS OKO ELH OOD) 
& 25 O = VIN] IPWIFVYER DOW ZVEIKO1D OWS = ies _ Bereta ees oa 
27 35 O = Yel) 


Oa aa RE UGN, 2 eS ee ee 





29 END 





SUBROUTINE SORTXYEX,Y,NPTS) 











1 — i a 
2 C Nee a Pi ng, at ge a EL, ge ye a ee. ee eS Ce as ans 
8 C THIS SUBROUTINE SORTS x INTO ASCEROING ORDER ee te 7 = 
a C 
———-—_3_____________DUEWSION xd 3NOrevesnoy at 
6 N = NPTS 
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22 10 CONTINUE I I I ene ges 
23 RETURN ———— Ka 
24 END 
w restates ciel 7 he he ae tee eet te ote mets saa 
® ea ee, . SoA at 4S ares ¥ Se Ses eee 2s 2 Sees ota ee, 





me 


a~ 















Se oe < ees tee re ae ES eee ee eee ee ee whe w FB e hae  & a8 

5 2 VCUZSXRI KAZ yXRH, VR1 yg VRZ2 eo YRHSCUTOF? o€RTOF 2.CUTOFH, ~ 
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cota ee Soe La V2" 1200) pV3XE 200) pV4X1 200) pV5X 12009 pBS¥I200), = tesa i Ectuts | 
8 2 V2V(200) oV3YC2C0) sV4Y¥ 0200) ,V5Yt200) eae oe ee Sa aS ON 

LLL MM ON sOOUTS NT 9NS1»NH pNP oT WoNX sKND pICOMP)K pNXMI] GNXHI2 oNXHI3, ae eee see ae 

& 10 1 NSTcoNST3eNST7 WNPFEISCDSTRAKE30) oM1 M7 eICOMPL,IPL. -~ 
1} 2 HL 
2 ( SASee TAS es ae ew en swinsas -Saceccaces so oecsascescoale cudccuclcecelel cL 

a 213s! CALL INTER(XONSYON oF e123} etn eee a i 
18 CALL INTERCKON, VI, IL gI2 hy CP 


15 CALL INTERCXON V2 91) oJ 2 ph gD) 
& 16 CALL INTERCKONSV3—5])2 512 ch ckD a 
1? CALL INTER(XON V4 J] IZ oA oF) 
18 IFITHUBSNE SO) CALL INTER(XON,V5,12,12 R06) 
@ 19 IF(IHUBeEQ.0) 62060 “ 
20 RETURN 
21 END 
— = ce ee 


@ i SUBROUTINE INTERS (Ap TD yT2 oC oD ol oF Hy P oR S2,7? “ 
2 COMMON /COORD/ XON 1700) »YONE 700) .XOF F1(2003,.VOFF¢20CDeSt700)~ 
3 i $1700), XTEST oXTE STI oXTESTZeVCL sYCU sVCL I> YCL2,YCUI, 
ae % 2 YCUZo KR1 oXR24 ARH, WR1g VR2 ep YRHGCUTOF } pCUTOF 2 eCUTOFH, sy 
® 7 


AR¢7C0) 

















6 COMMON /SOLUT/ VBARE700),VBAROI2CO),VINCE 700) » ¥XINCI200) w 
@ 7 1 VY IAC (200) gRHOBI7(C) eRBORT C700) RHOPO(2003, _o© : 
8 2 VCORE700},RBOOTI2ZCOD, VRE C200) pVRECOMEZ00)« 
sel 2@e 2 Siete VXC CH 200) pVYCOMI 200) o THETAL2003 ePSCPICITO0? Ste ae ee 
e 10 ® PSOFT (700) sCMACHE 700) »XMACHE700 DCP 10700) eCPCET00). a 






il RHOJt7O 


eo az DIMENSION At700) a 
14% CALL INTERGA,YON, 11,12,C,0) 
15 CALL INTERIM oVCOM, T1,I2,C,£) ge oe Soe cee! 
e@ 16 CALL INTERGA oVINC, Tl el2eC oF) ot 
17? CALL INTERGCA SRBORI al by 125C pH?) 
18 CALL INTERGA oVBAR, T1,12.C,P) 
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Figure 5. - Inlet and nozzle geometry 
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(a) Single passage nozzle 
At any point 
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1 2 of the three control stations 
Veo free stream velocity 
a direction of free stream velocity relative 
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(b) Two passage nozzle 


Figure 7 - Combined solution, V 
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Figure 8 - Graphic output of SCIRCL 
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Figure 1] - Vector flow field from nozzle 
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calculated. 
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Exponents P and Q calculated 


(e) Optional point (Xz Y3) and 
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P and Q calculated. 





Figure 12. - Sketches for SCIRCL input. Bisuperellipse options 
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(g) Curvature at either point 2 or 
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tated curvature at endpoint 
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Figure 12. - Concluded 
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Figure 13 - Sketches for SCIRCL tnput options except bisuperellipse 
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(a) Title, control, and rake data cards. 


Figure 14. - SCIRCL input layout. 
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(a) Title, control, and rake data cards. 


Figure 14, - SCIRCL input layout. > 
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Figure 15. - NOZZLEC input form 










MAIN 
SCIRCL 
4B,W,B,R 
5:R(P.S, SC, L) 


PLOXIS 
(P,S,N) 


ae STRAIT 


ale 


(P, L, SC) 
FRSTSH 
-{ so -_ 
ECHO f+ ERTRAN | TEST 
PRELPS 
DRAW 
(S,L) 
ELIPSE 
17:W 
17:W 


SUPERC 


WPUNCH 
(L) WRIXY 


lee) (ELE) BEE 


CUBIC 


MIRROR 


5 - 1/0 Unit No 
B - Rewind 
W - Write 
R - Read 
(a) SCIRCL 


Figure 16. - Call Sequences. 


FNTR PA DNTRPC 
LENTRP_ | 


XYCALC 


SGEN | 
DZTRP 


FNTR PC 


Calcomp Routines Referenced 


(L) - Line 

(P) - Plot 

(S) - Symbol 
(SC) - Scale 

(N) - Number 


5,6,7: RW 
MAINT 


(b) 24Y 
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12 13, 14B, R, W 
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SORTXY 


(c) NOZZLEC. 
Figure 16. - Concluded. 
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